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ABSTRACT 


For  cylindrical  shells  non-axi symmetrical  loading  is 
generally  considered  in  the  form  of  a  periodic  Fourier  series. 
Since  each  term  of  the  series  must  satisfy  the  boundary  con¬ 
ditions,  the  solution  is  simplified  by  the  use  of  stiffness 
coefficients.  To  permit  a  method  of  solution  closely  following 
the  classical  stiffness  matrix  method  expressions  for  stiffness 
influence  coefficients  for  all  the  harmonics  of  loading  greater 
than  zero  are  presented  as  general  algebraic  expressions  in 
closed  form.  Two  sets  of  expressions  were  obtained  by  introducing 
approximations  into  the  classical  theories.  The  range  of 
applicability  of  these  sets  is  given  based  on  criteria  involving 
the  shell  parameters  and  the  harmonic  number  of  the  edge  load¬ 
ing  considered.  A  further  simplification  is  possible  for 
evaluating  membrane  stresses  for  very  large  values  of  the  load¬ 
ing  harmonic. 

The  values  of  coefficients  and  the  stress  resultants 
obtained  for  two  examples  are  compared  with  those  obtained 
using  Flugge's  equations  for  these  shells.  It  is  concluded  that 
the  use  of  these  coefficients  constitute  a  practical  procedure 
for  the  solution  of  these  types  of  problems. 
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CHAPTER  I 


INTRODUCTION 

1 . 1  PURPOSE  OF  THE  STUDY 

An  unsymmetri cal  state  of  stress  is  obtained  in  shells 
of  revolution  such  as  water  tanks,  nuclear  reactors,  pressure 
vessels,  silos  and  chimneys  by  one  of  these  causes: 

1.  arbitrary  loading  varying  along  the  directrix  as 
in  the  case  of  wind  and  seismic  actions, 

2.  the  supporting  reactions  being  ' unsymmetri cal '  or 
discontinuous  as  in  the  case  of  a  shell  supported  on  isolated 
col umns , 

3.  the  middle  surface  being  incomplete. 

The  symmetrical  state  of  stress  results  from  particular  cases 
of  these  conditions. 

The  classical  method  of  analysis  for  arbitrary  loading 
is  to  express  the  loading  as  a  Fourier  series,  analyse  for  each 
harmonic  separately  and  then  superimpose  the  effects.  The  number 
of  terms  considered  in  this  series  must  be  sufficient  to  give  the 
desired  degree  of  convergence.  An  example  of  this  procedure  is  the 
consideration  of  wind  and  seismic  loading.  With  regard  to  wind 
action  the  normal  component  pw  of  the  equivalent  static  load  is 
expressed  as 

pw  =  9w  k(4,) 


1 


2 


where  k(<|>)  is  the  shape  coefficient  and  gw  is  the  basic  dynamic 
pressure  which  depends  on  factors  affecting  wind  pressure  distribu¬ 
tion.  From  studies  on  cylindrical  models  in  wind  tunnels  several 
investigators  have  expressed  the  law  of  variation  kfc  )  for  smooth 
and  rough  surfaces  as  a  trigonometric  series  in  the  form, 

k(<J>)  =  aQ  +  a-jcosf  +  a2Cos2t  +  a3cos3<f>  +  a^cos4f  +  ... 

The  usual  procedure  in  design  is  to  consider  only  the  first 
harmonic  a^cosf  ignoring  the  higher  harmonics.  The  analysis 
is  accomplished  by  the  beam  theory  treating  the  structure  as 
a  cantilever  beam  or  by  the  membrane  theory  of  shells.  This 
procedure  is  satisfactory  only  if  the  edges  fully  satisfy  the 
required  boundary  conditions.  According  to  Hruban  (1)  for 
certain  stress  resultants  such  as  the  bending  moments,  the 
higher  harmonics,  especially  the  second  one,  may  represent 
very  significant  components. 

An  example  of  higher  harmonic  loading  due  to  intermittent 
supporting  reactions  is  the  case  of  a  cooling  tower  supported 
on  eight  columns.  Here  the  reactions  when  expressed  as  a 
Fourier  series  will  contain  the  harmonics  8,  16  ,  24  etc. 

With  the  advent  of  the  digital  computer,  structural 
analysis  may  be  conveniently  done  by  matrix  methods  using  in¬ 
fluence  coefficients.  In  the  literature,  the  analysis  for 
symmetrical  loading  is  usually  performed  using  influence  co¬ 
efficients  which  have  been  given  in  an  explicit  manner.  For 
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unsymmetri cal  bending  there  are  no  correspond! ng  explicit 
expressions.  Although  one  can  theoretically  write  the  exact 
expressions  based  on  a  rigorous  theory,  this  procedure  involves 
the  solution  of  eight  simultaneous  equations,  expressed  in 
general  terms  of  geometry  of  the  shell  surface  and  the  physical 
constants.  Moreover,  the  final  expressions  for  influence  co¬ 
efficients,  when  written  in  an  explicit  manner,  will  be  too  unwieldy 
to  be  of  practical  use.  For  these  reasons  the  expressions  for 
influence  coefficients  have  to  be  established  based  on  approximations. 
Therefore  the  approximate  methods  of  analysis  proposed  in  the 
literature  are  reviewed  in  detail  in  the  next  section. 

The  purpose  of  this  study  is  to  investigate  the  possibility 
of  establishing  these  influence  coefficients  in  closed  form  for 
practical  use. 

1 .2  A  BRIEF  HISTORICAL  REVIEW  OF  SHELL  THEORY 

The  general  theory  of  thin  shells  was  first  initiated 
by  Love  (2)  in  1890.  Flugge  (3)  in  1932  presented  the  basic 
equations  of  the  theory  in  the  form  mostly  used  today.  Because 
of  the  complicated  nature  of  the  rigorous  theories  proposed  by 
Love  and  Fliigge,  simplifications  were  proposed  by  Finsterwal der  (.4) 
in  1932,  Donnell  (5)  in  1934  and  Shorer  (6)  in  1935.  It  was  not 
until  after  1940  that  Donnell's  theory  gained  wide  recognition. 

In  establishing  his  equations  Donnell  rigorously  discarded  the 
lesser  important  terms.  The  Donnell  theory  was  apparently  derived 
independently  by  Von  KaVman  (7)  in  1941,  Jenkins  (8)  in  1947  and 
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Vlassov  (9)  in  1947. 

The  problem  of  unsymmetrical  bending  of  closed  cylinders 
was  first  investigated  by  Schwerin  (10) in  1922  followed  by 
Miesel  (11 ) i n  1930.  Hoff  (12)  in  1954  gave  an  explicit  solution 
for  Donnell's  equations.  In  1957  Holand  (13)  proposed  explicit 
solutions  using  a  modified  form  of  Flugge1 s  equations.  An 
approximate  method  of  analysis  with  possible  simplifications 
was  given  by  Flugge  (j4)  in  1960.  Goldenveiser  (15)  in  1961 
made  a  detailed  study  of  the  mathematical  aspects  of  the  shell 
theory  and  proposed  simplified  solutions. 

1 .3  REVIEW  OF  EXISTING  METHODS  OF  ANALYSIS 

A  detailed  examination  of  the  more  important  works 
considering  unsymmetrical  bending  of  closed  circular  cylindrical 
shells  follows.  Goldenveiser  and  Flugge  have  presented  discussions 
on  the  nature  of  the  shell  characteristic  equations,  while  Hoff 
and  Holand  have  proposed  solutions  to  these  equations. 

GOLDENVEISER 

2  2  2 

Making  the  assumption  that  a<<l  where  a  =  h  /3r  , 
Goldenveiser  obtains  the  following  eighth  order  characteristic 
equation. 

x8  -  4m2 *x6  +  [6m4  -  (8-2u2)m2  +  (l-u2)]x4 

2 

-  4m2 (m2  -  if  x2  +  m4  (m2  -  l)2  =  o  (1-1) 

For  values  of  m  ?  1  ,  this  equation  reduces  to  the  following 


forms . 
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For  m  =  1 


(1  -  u2)a4  =  0 

A4  +  (1  -o2)  =  0 


(l-2a,b) 


For  1  <  m  <  a  2 


ft  2  \  4  ,  28 

( 1  -  u  )  A  +  a  m  =  o 


x4  +  (1  -  u2)  =  o 


(1 -3a  ,b) 


2 


a 


For  m  ^  a  2 


(A2  -  m2)4  +  (1  -  u2)  A4  =  o 


(1-4) 


2 


a 


For  m  a  z 


(a2  -  m2)4  =  o 


(1-5) 


Depending  on  the  value  of  m  each  form  of  these  eqations  will 
have  eight  complex  roots,  A^  to  Ag  ?  expressed  in  terms  of  real  and 

imaginary  parts,  XlsX2  and  ^ ,y  , respectively .  EQUATIONS  (l-2a)  and 
( 1  - 3a )  are  used  to  obtain  the  four  roots  expressed  in  terms  of 
Xg,  yg  which  are  small  in  magnitude  and  EQUATIONS  ( 1  - 2b )  and 
( 1  - 3b )  are  used  to  obtain  the  four  roots  expressed  in  terms  of 
Xy'  y-|  which  are  large  in  magnitude.  EQUATION  (1-4)  as  well 
as  EQUATION  (1-5)  gives  all  the  eight  roots  expressed  in  terms 
of  x-j  >  y-j  and  x2>  y2>  xi  an^  x2  *Dein9  same  order  of  magnitude. 

The  stresses  and  displacements  at  any  point  in  the  shell 
are  expressed  in  terms  of  a  potential  function  $  as 


* 


' 
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For  l^m<a2the  X\~  solutions  are  termed  the  edge 
effects'  because  of  their  highly  damped  nature  and  the 
solutions  as  the  'basic'  because  of  their  lightly  damped  nature. 
Form^a  &  this  classification  does  not  apply  since  the  rate  of 
damping  is  similar  in  both  the  cases. 

To  consider  the  boundary  conditions  for  any  particular 
problem,  it  is  necessary  to  evaluate  the  eight  constants  of  in¬ 
tegration  C-|  to  Cg.  Goldenveiser  proposes  two  methods  for 
evaluating  these  constants,  one  applicable  for  nKa  Jand  the  other 
for  m?a  2' 

For  m<a  2  the  stress  system  corresponding  to  the  basic 
state  and  the  edge  effects  can  be  determined  separately  using 
EQUATIONS  (1-2)  or  (1-3)  and  then  superimposed.  According  to 
this  any  tangential  displacement  or  force  T  can  be  expressed  as, 

T  =  TbaS(y,C1,C2,C3.C4)  +  Ted9(y,C5,C6,C7,C8)  (1-7) 

In  a  similar  manner  any  non  tangential  displacement  or  force 
can  be  written  as, 

S  =  SbaS(y1,CrC2,C31C4)  +  Sedg(y,C5>C6)C7,C8)  (1-8) 

In  the  basic  state  of  stress,  tangential  forces  are  large  in 
comparison  with  transverse  forces  and  moments.  In  the  edge  effect, 
the  strain  in  the  middle  surface  of  the  shell  is  determined 
primarily  by  the  normal  deflection  and  rotation.  It  therefore 
follows  that  on  the  right  hand  side  of  EQATION  (1-7)  the  second 
term  is  much  smaller  compared  to  the  first  term.  Neglecting 


' 
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this  term,  EQUATION  (1-7)  can  be  written  as, 


T  =  TbaS(y,C, ,C9>C,,C.) 


1  2 ’3’  4 


(1-9) 


Using  this  equation  the  constants  C-j  to  can  be  determined  from 
a  set  of  four  tangential  boundary  conditions.  Recognising  the 
highly  damped  nature  of  the  x-|  -  solutions,  EQUATION  (1-8)  can 
be  written  for  the  top  and  bottom  edges  separately  as, 


As  C-j  to  have  been  already  found,  the  remaining  constants  to 
Cg  can  be  determined  using  EQUATIONS  (1-1 Oa , b )  from  a  set  of 
four  non  tangential  boundary  conditions. 


t 


Form^aethe  problem  is  governed  by  the  single  eighth 
order  EQUATION  (1-4)  or  (1-5).  Since  both  x-|  and  X2“s°l utions 
are  highly  damped  the  boundary  conditions  may  be  applied  in¬ 
dependently  of  each  other  at  the  two  edges.  The  potential  function 
$  is  expressed  as 


$  =  |e"xl  [C-jCosy-jy  +  CgSiny-jy]  + 
e"X2  [C5cosp2y  +  CgSinp2y]J  cos  m* 


(1-11) 


There  are  now  only  four  constants,  C -j C ^ , C ^  and  Cg,  to  be  determined 
from  a  set  of  four  boundary  conditions  at  the  edge  where  the 
loads  are  applied. 
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For  m»a  a  that  is  for  very  large  values  of  m  the 
stress  system  can  be  subdivided  into  generalised  plane  state 
of  stress  and  plate  bending.  The  tangential  displacements 
and  forces  are  given  by  the  equations  of  the  generalised 
plane  state  of  stress.  The  non  tangential  displacements 
and  forces  are  given  by  the  equation  of  plate  bending.  The 
equation  of  equilibrium  for  the  shell  element  are  given  by  these 
two  groups  of  uncoupled  equations  each  leading  to  a  single 
characteristic  equation  of  the  following  form. 


0~12) 


o 


Thus  for  very  large  values  of  m  the  shell  element  behaves  as 
if  it  were  developed  into  a  thin  plate. 


FLUGGE 


The  eighth  order  character!- Stic  equation  given  by  Flugge  is 


8  0/02  *6  ,  T  !  2  ,  2  /  2  -i\,4 

X  -  2  (2m  -  u)x  +  1  -  u  +  6m  (m  -  1 )  X 

-  k 


+ 


o 


(1-13) 


For  small  values  of  m  this  equation  can  be  split  into 


two  independent  fourth  order  equations  as, 
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A4  -  2(2m2  -  u)a2  + 


i  2  ,  c  2/  2 

-  u  +  6m  (m  -  1 ) 


k 


=  o 


1  -  u2  +  6m2 (m2  -  1  ) 

L  k 


A  9  A  ? 

A  -  2m  [2ni  -  (4  -  u)nr  + 


(2  -  u)]  A2  +  m4 (m2  -  1)  =o  (1 -14a ,b) 

EQUATION  (1-1 4a )  gives  the  larger  roots  expressed  in 
terms  of  x-|  >  y-i  and  EQUATION  (1-1 4b )  gives  the  smaller  roots 
expressed  in  terms  of  x£>  the  roots  correspond!* ng  closely 
with  those  obtained  from  EQUATION  (1-13).  The  large  roots  are 

From 


_1  /4  1 /4 

of  the  order  k  and  the  small  roots  of  the  order  k 


the  table  of  characteristic  coefficients  for  the  displacements 
and  the  stress  resultants  (TABLE  1,  p.  230,  Tef.  14)  the  states 
of  stress  corresponding  to  EQUATIONS  (1-1 4a , b )  can  be  determined 
separately.  In  this  table,  the  quantities  a-  are  of  the  orders 

1 /4  _1 /4 

k  and  k  for  the  large  and  small  roots  respectively  and 


1 


/2 


3-  of  the  orders  k  and  1.  It  is  seen  that  for  the  state  of 

vJ 

stress  and  displacements  given  by  EQUATION  (1-1 4a )  the 

membrane  quantities  are  very  small  compared  to  the  bending 

quantities  and  vice  versa  for  the  state  given  by  EQUATION  (1-1 4b ) . 

This  means  that  the  state  of  stress  corresponding  to  EQUATION 

(1-13)  can  be  found  approximately  by  determining  separately  the 

two  states  of  stress  as  mentioned  above.  Therefore  EQUATION  (1-1 4a ) 

satisfy  a  pair  of  non  tangential  boundary  conditions  Mx  or  W 

or  SY  or  w  and  EQUATION  (1-1 4b )  will  satisfy  a  pair  of  tangential 

boundary  conditions  N  or  u  and  T  or  v.  The  boundary  conditions 

x  x 
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corresponding  to  EQUATION  (1-1 4b )  are  those  considered  in  the 
membrane  analysis.  Therefore  this  part  of  the  solution  is  only 
an  improved  form  of  the  membrane  solution  and  the  membrane  sol¬ 
ution  can  be  used  as  an  approximation  for  the  solution  of 
EQUATION  (1 -14b) . 

For  large  values  of  in,  the  splitting  of  EQUATION  (1-13) 
in  the  above  manner  is  not  valid  and  values  of  x-|  approach  x2 
and  values  of  y-j  approach  Because  of  the  highly  damped 
nature  of  both  x-j  and  x2  solutions  the  effect  of  far  edge 
can  be  neglected,  thus  requiring  the  consideration  of  only 
four  boundary  conditions. 

HOFF 

Donnell  ' s  characteristic  equation  is  expressed  as 

(a2  -  m2)4  +  1  A4  =  o  (1-15) 

k' 

Using  this  equation  Hoff  (12)  obtained  expressions  for  the  middle 
surface  displacements  in  terms  of  eight  constants  which  are 
evaluated  from  the  boundary  conditions.  This  was  accomplished 
by  finding  solutions  for  the  root  of  the  characteristic  equation 
and  the  relations  among  the  middle  surface  displacements. 

HOLAND 

In  an  attempt  to  obtain  explicit  solutions  to  the  roots 
of  Flugge's  characteristic  equation,  Holand  modified  his  equation 
in  the  following  form, 


11 


■i 


4m2 (in2  -  1)  m/m2  -  lA2  +  m4(m2  -  l)2 


0 


(1-16) 


The  equation  can  also  be  written  in  the  form 


(1-17) 


The  value  of  the  roots  as  given  i n (SECTION  A.l  of  the  Appendix) 


depend  on  the  parameter  g  which  is  defined  as  4 fk'mjrn^  -  1 
For  m  =  o  or  1 , 


and  X2  ~  ^2  =  0 


X1  =  v\  = 


The  roots  are  therefore  the  same  for  both  m  =  o  and  m  =  1 . 

For  small  and  large  values  of  g,  Holand  gives  approxim¬ 
ate  relations  for  the  roots.  He  further  suggests  corrections 
to  the  roots  obtained  from  EQUATION  (1-17).  But  for  all 
practical  purposes  where  the  shell  is  thin,  the  approximate 
solutions  are  very  close  to  the  exact  ones. 

Holand  presents  a  discussion  of  the  accuracy  of  the 
solution  of  the  different  characteristic  equations,  proposed 
by  many  investigators.  If  (m  -  1)  is  replaced  by  m  in 
EQUATION  (1-17) ,  Donnell 's  equation  is  obtained.  Neglecting 
terms  which  are  of  no  importance  to  the  design  of  thin  shells, 
he  obtains  simple  expressions  for  his  characteristic  coefficients. 


Explicit  expressions  for  influence  coefficients  neglecting 


the  effect  of  far  edge  are  also  obtained  by  Holand  .  He 


- 
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considers  a  set  of  mixed  boundary  conditions  expressed  by  the 

group  of  edge  quantities  M  N  and  v  and  w.  For  m  =  1 ,  he  obtains 

x’  x’ 

the  solutions  in  two  parts,  a  polynomial  solution  and  a 
damped  exponential  solution.  He  concludes  that  the  membrane 
theory  can  be  used  in  place  of  the  polynomial  solution. 

Explicit  expressions  for  both  stresses  and  displacements  are 
given  for  each  part  of  the  solution. 

1 .4  AUTHOR’S  CLOSURE  TO  LITERATURE  REVIEW 

The  differential  equations  obtained  in  the  general  theory 
of  shells  as  encountered  in  the  literature  are  quite  complicated 
even  when  based  on  a  consistent  first  approximation  as  in  the 
sense  conceived  by  Love  (2).  To  obtain  solutions  it  has  been 
found  necessary  therefore  to  introduce  simplifications.  The 
review  of  literature  points  out  to  two  types  of  simplifications. 

In  the  first  type,  approximate  solutions  result  from  a  con¬ 
sideration  of  the  states  of  stress.  It  is  then  necessary  to 
proceed  to  a  classification  of  the  possible  states  of  stress 
and  for  each  state  of  stress,  consider  which  simplifications 
or  approximations  are  possible.  This  does  of  course  result 
in  a  loss  of  the  general  character  of  the  theory.  The  approximate 
methods  of  Goldenveiser  and  Flugge  belong  to  this  category. 

The  second  method  is  based  on  solutions  to  approximate 
eighth  order  characteristi c  eqations.  Hoff's  method  based  on 
Donnell's  equation  and  Hoi  and' s  simplification  of  Flugge' s 


'J'j? 


13 


equation  belong  to  the  second  category. 

Only  Goldenveiser  (15)  has  done  a  detailed  study  of 
cylindrical  shells  with  bending  resistance  under  arbitrary 
loads.  Within  the  mathematical  framework  he  defines  the 
region  of  applicability  of  his  solutions  with  an  estimate 
of  errors  in  the  roots.  He  has  not  made  any  quantitative 
study  regarding  the  error  in  the  final  stress  resultants. 

Flugge  (14)  has  suggested  similar  simplifications 
but  has  not  presented  any  specific  criteria  to  decide  when 
the  problems  can  be  split. 

By  introducing  various  assumptions  different  investigators 

have  obtained  simplified  characteristi c  equations.  TABLE  1.1 

gives  a  comparison  of  the  different  characteristic  equations. 

Flugge' s  equation  is  the  most  accurate  one  based  on  a  technical 

theory  of  shells  in  that  he  neglects  only  the  fifth  and  higher 

powers  of  h_  in  the  development  of  his  theory.  The  other 
r 

characteristic  equations  are  therefore  compared  with  that 
of  Flugge. 

The  difference  between  the  equations  of  Holand  and 

6  2 

Flugge  appears  in  the  coefficients  of  A  and  A  .  For  the 

usual  type  of  shells  k  is  of  the  order  10“4  to  lCf7.  As  such 

4  6  2 

the  coefficient  of  A  is  very  large.  The  terms  A  and  A  are 

of  little  importance.  Holands'  equation  is  therefore  the  best 

approximation  given  so  far  to  Fliigge's  exact  equation. 


14 

Morley  (16)  makes  assumptions  which  are  tantamount  to 

p 

taking  o  =  1  in  Flugge's  equation  ,  except  in  the  term  (1  -  u  )/k. 
This  indicates  that  the  roots  are  insensitive  to  variation  in  o. 
Morley's  computed  values  of  A  are  only  slightly  in  error. 

Bijlaard's  (17)  equation  is  accurate  for  large  values  of  tn. 

Donnell's  equation  is  widely  used  because  of  its  simplicity. 

In  this  equation  only  the  higher  powers  of  m  are  retained.  Hoff  (18), 

Morley  (16)  and  Kempner  (19)  conclude  that  the  accuracy  of  Donnell's 

equation  improves  as  m  increases.  However,  an  investigation  by 

Col  bourne  (22)  shows  that  when  rn  increases,  the  real  parts  of  the 

roots,  anc^  X2 9  approach  the  exact  values  but  the  imaginary  parts, 

and  , do  not  agree  with  the  exact  values.  Thus  the  effect  of  the 

error  in  the  quantities. P2  and  p2,on  final  stress  resultants, 

for  high  values  of  m,  needs  investigation  and  is  discussed  in 

CHAPTER  V.  Hoff  points  out  that  Donnell's  equation  should  not  be 

used  when  m  is  less  than  4.  Biezeno  and  Grammal  (20)  make  a  further 

2  2 

approximation  to  Donnell's  equation  by  neglecting  — 

compared  to  6m^.  The  solution  is  then  independent  of  k. 

It  is  seen  that  the  values  of  the  roots  in  the  different 

2 

characteristic  equations  depend  on  'm'  and  1  -  u  »  According 

k 

to  Hoff  (16)  when  1  -  u2  <2500  the  shell  is  relatively  thick 

k 

and  the  simplifying  assumptions  made  in  deriving  the  basic 

2  7 

equations  may  invalidate  the  basic  theory.  When  1  -  u  >2.5x10 

k 

it  represents  the  thinnest  practical  shell.  r 
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A  different  kind  of  inaccuracy  results  due  to  an 
inconsistency  in  the  theory  when  one  uses  approximate  equations 
like  that  of  Donnell  or  Bijlaard.  For  m  =  1  these  equations 
do  not  give  the  four  zero  roots.  These  are  essential  for  rigid 
body  displacements  required  from  equilibrium  considerations.  Any 
characteri Stic  equations  will  have  the  necessary  number  of  zero 
roots  only  if  the  elasticity  relations  are  chosen  in  the 
proper  way.  By  changing  the  relations  even  within  the  frame¬ 
work  of  the  exact  theory  of  shells  one  can  lose  the  zero 
roots  and  obtain  large  errors  in  the  final  results. 

On  reviewing  the  literature  it  is  seen  that: 

1.  approximate  solutions  even  where  they  exist  have  not 
been  thoroughly  investigated  with  regard  to  their  effects  on 
the  final  solutions  and 

2.  no  explicit  expressions  exist  for  stresses  and  dis¬ 
placements  leading  to  the  development  of  stiffness  influence 
coefficients. 

Therefore,  the  objectives  of  the  present  study  can  be  stated 
as  follows: 

1.  to  study  the  effect  of  variation  in  the  shell 
parameters  and  the  harmonic  number  of  the  edge  loading  on  the 
behaviour  of  the  shell, 

2.  to  arrive  at  values  of  stiffness  influence  coefficients 
for  particular  values  of  the  shell  parameters  and  harmonic  number 
of  the  edge  loading,  using  the  exact  theory, 
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3.  to  establish  approximate  closed  form  expressions  for 
stiffness  influence  coefficients  for  all  harmonics  of  loading 
and  to  establish  the  range  of  applicability  of  these  expressions, 

4.  to  solve  certain  specific  problems  using  these 
expressions  and  to  compare  the  results  with  those  obtained  by 
the  exact  theory. 


TABLE  1.1 


COMPARISON  OF  CHARACTERISTIC  EQUATIONS  FOR 
FOURIER  SERIES  REPRESENTATION  CIRCUMFERENTIALLY 


Flugge 

A8  +  2(2m2-u)A6  +  [1-u2  +  6m2 (m2-l )  A4 

k 

-2m2[2m4  -  (4-u)m2  +  (2-u)]a2  +  m4(m2--l)2 

=  0 

Holand 

A8  -  4m  / m 2 - 1  A8  +  [1-u2  +  6m2(m2-l)]A4 

u  k 

-4m2 (m2-l )m J m2-l  A2  +  m4(m2-l)2 

=  0 

Morley 

A8  -  2  (2m2-l)A6  +  [1-u2  +  6m2(m2-l)]A4 

k 

-2m2(2m4-3m2+l )A2  +  m4(m2-l)2 

=  0 

Bi  jlaard 

A8  -  2  (2m2)A6  +  [1-u"1  +  6m2(m2-l )  -  u(l-u 

k 

)nP]A4 

-2m2  [2m4-^/2X7+u)m2]A2  +  m4  (m4-2m2) 

=  0 

Donnel 1 

A8  -  2(2m2)A6  +  [1-u2  +  6m2(m2)]A4 

k 

-  2m2(2m4)A2  +  m8 

=  0 

Donnel 1 
modif i ed 
by  Biezeno 
&  Gramme! 

A8  -  4m2A6  +  6m4A4  -  4m6A2  +  m8 

=  0 

CHAPTER  II 


INFLUENCE  COEFFICIENT  PROCEDURE 

2.1  APPLICATION  OF  MATRIX  METHODS  TO  STRUCTURAL  PROBLEMS 

Structural  problems  which  usually  consist  of  an  assemblage 
of  many  parts' tend  to  be  complex  in  nature.  The  true  structure 
must  generally  be  replaced  by  an  idealised  approximation  or 
model  suitable  for  mathematical  analysis.  Matrix  methods  using 
influence  coefficients  permit  a  more  detailed  idealisation  of 
structures  than  was  possible  by  the  earlier  classical  procedures. 

Matrix  methods  can  be  classified  as  force  methods  using 
flexibility  influence  coefficients  and  displacement  methods 
using  stiffness  influence  coefficients.  Such  influence  co¬ 
efficients  relating  forces  and  displacements  reduce  the  given 
problem  to  a  mathematical  formulation  in  terms  of  a  set  of 
linear  algebraic  equations  rather  than  differential  equations. 

The  solution  of  a  large  number  of  simultaneous  equations  is  no 
more  a  problem  with  the  advent  of  high  speed  electronic  computers. 
Matrix  methods  have  therefore  come  to  be  well  recognised  in 
structural  analysis.  Both  force  and  displacement  methods  can 
be  applied  to  a  variety  of  complicated  engineering  problems. 

The  following  basic  principles  are  applied  while  setting 
up  both  types  of  influence  coefficients. 


18 


19 

1.  Equilibrium  of  forces 

2.  Compatibility  of  deformations 

3.  Constitutive  equations 

2 . 2  FLEXIBILITY  AND  STIFFNESS  INFLUENCE  COEFFICIENTS 

Assume  an  arbitrary  linearly  elastic  structure  supported 
against  rigid  body  motion  and  subjected  to  forces  . .  •  • 

acting  at  nodes  1,  2,  . ...,  n.  Let  the  correspondi ng  displace¬ 
ments  be  d^  ,  d^, . ,  d  .  The  terms  force  and  displacement 

are  used  here  in  a  general  sense.  In  matrix  notation  the 
relation  between  the  forces  and  the  displacements  can  be  expressed 
as , 

fd]  =  [C]  [X]  (2-1) 

where  d  is  the  column  vector  {d-j,  d^,....,  d  }  and  [X]  is 
another  column  vector  {x-|>  ^ 

The  square  matrix  [C]  is  known  as  the  matrix  of  flexibility 

influence  coefficients.  The  element  c. .  is  defined  as  the  dis- 

ij 

placement  at  the  node  i  due  to  unit  force  applied  at  the  node  j  • 
EQUATION  (2-1)  represents  load-deflection  equations  for  the 
structure.  The  forces  in  terms  of  deflections  can  be  solved 
from  these  equations. 

For  the  same  problem,  the  relations  given  by  EQUATIONS  (2-1) 
can  also  be  restated  as 


[x]  =  [K]  [d] 


(2-2) 
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where  [K]  is  known  as  the  matrix  of  stiffness  influence  co¬ 
efficients.  The  element  k..  is  defined  as  the  force  at  the 

*  vJ 

node  i  due  to  unit  displacement  at  the  node  j.  EQUATIONS  (2-2) 
represent  the  equilibrium  conditions  for  the  structure.  The 
displacements  in  terms  of  forces  can  be  solved  from  these 
equations . 

Because  of  the  reciprocal  theorem  both  [C]  and  [K]  are 
symmetric. 

2-3  DIRECT  STIFFNESS  METHOD 

While  establishing  the  flexibility  matrix  the  structure 
is  assumed  to  be  supported  against  rigid  body  motion,  a  condition 
not  necessary  for  the  stiffness  matrix.  For  the  latter  case 
the  structure  can  be  free  to  move  as  a  rigid  body  when  the  edge 
displacement  is  applied.  The  stiffness  matrix  [K]  thus  obtained 
is  called  the  'direct  stiffness  matrix.' 

The  symmetric  direct  stiffness  matrix  [K]  is  singular. 

The  elements  in  any  column  represent  a  stress  system  in  equilibrium. 
The  singular  nature  of  [K]  can  be  removed  by  imposing  boundary 
conditions  on  the  problems.  The  advantage  of  the  direct  stiff¬ 
ness  matrix  method  is  that  for  a  whole  assembly  of  structural 
elements,  the  total  structural  stiffness  matrix  can  be  easily 
formed  by  superposition  of  the  individual  stiffness  matrixes, 
irrespective  of  the  boundary  conditions.  The  boundary  con¬ 
ditions  are  considered  only  in  the  actual  solution  of  the 
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system  of  equations.  This  permits  the  consideration  of  different 
boundary  conditions  while  the  total  structural  stiffness  matrix 
remains  unaltered. 


2-4 


SOLUTION  PROCEDURE  FOR  DIRECT  STIFFNESS  METHOD 


The  total  structural  stiffness  equations  are  represented 


by 


(2-3) 


[Xj  =  [K]  [dj 

where  [x]  is  the  total  set  of  nodal  forces  and  [d  \  ,  the 
corresponding  total  set  of  nodal  displacements.  Let  [d  j  represent 
the  unknown  nodal  displacements  and  [d  “1  the  known  nodal  displacements 
due  to  the  boundary  conditions.  Correspondingly  [x  ]  can  be  separated  as 


the  applied  loads  [X^|  and  the  reactions  [X^]  .  Thus, 


Pal 

tdbl 


Tg 

9* 

K 

K  ' 

i 

aa 

ab 

]xbi_ 

/ba 

Kbb 

I— 

(2-4) 


iXal 


K. 


aa 


[dal  +  [Kab]  fdbi 


K 


ba 


K1 


+  K 


bb 


M 


(2-5) 

(2-6) 


From  EQUATION  (2-5) 

[  dal  =  V'  [KJ  -  ^  [dbf  • 

If  fdb]is  a  null  matrix  correspondi ng  to  a  set  of  boundary  dis¬ 
placements  specified  as  zero  ,  then 
EQUATION  (2-7)  becomes, 


(2-7) 


-1 


fdal  =  tKaa3  tXal 


(2-8) 
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Values  of  the  elements  in[X  ]  can  be  found  by  substituting  the 
values  of  the  displacements  Pd  1  obtained  from  the  EQUATIONS  (2-7) 
into  EQUATION  (2-6).  If  [d^]  is  a  null  matrix, 


P 


b .. 


-1 


[K.  ]  [K  ]  f  X  ] 

L  Ua  1  L  PtPi  1  L  a  -I 


(2-9) 


baJ  L  aaJ  a; 

When  all  the  nodal  displacements  are  known,  each  structural 
element  can  be  assumed  to  be  subjected  to  an  equivalent  system 
of  nodal  forces  designated  byp(  j,  where 

f>y  =  [K]  [d]  ,  (2-10) 

M  is  the  member  stiffness  matrix  and  [d]  the  actual 
edge  displacements  for  the  members.  The  internal  forces  are 
related  to [xj by  a  relation  directly  obtained  from  the 
derivation  of  the  member  stiffness  matrix. 


2-5  APPLICATION  OF  DIRECT  STIFFNESS  METHOD  TO  SHELL  PROBLEMS 

The  procedure  that  is  used  in  the  classical  shell  analysis 
is  to  consider  the  general  solution  which  consists  of  a  homo¬ 
geneous  part  and  a  particular  part.  The  general  solution  must 
satisfy  the  specified  boundary  conditions.  Using  the  stiffness 
matrix  given  by  EQUATIONS  (2-3)  the  homogeneous  part  of  the 
solution  can  be  written  in  terms  of  the  unknown  ^ 
edge  displacements.  The  particular  solution  or  the  membrane 
solution  is  superimposed  on  the  homogeneous  solution.  The 
resulting  general  solution  is  made  to  satisfy  the  specified 
boundary  conditions  and  the  unknown  displacements  are  evaluated. 


, 
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As  explained  in  the  previous  section,  the  internal  forces 
corresponding  to  the  homogeneous  solution  can  b 


can  be 


determined.  To  get  the  net  values  of  the  displacements  and 
the  stress  resultants  the  particular  solution  is  to  be 
superimposed  on  the  homogeneous  solution.  The  classical 
shell  analysis  therefore  does  not  give  the  net  values  of  the 
edge  displacements  directly. 

An  alternative  method  which  is  analogous  to  the  slope 
deflection  procedure  for  beam  elements  permits  a  matrix  formulation 
to  the  shell  problem.  In  the  local  co-ordinate  system  in  which 
the  final  displacements  are  directly  obtained  the  total 
structural  stiffness  equations  for  a  single  shell  with  no  surface 
loading  is  given  by 


(2-11) 


In  a  complete  analysis  the  surface  loading  also  has  to  be 
considered.  Let  £X  !]  and  j^d'^^be  the  edge  forces  and  dis¬ 
placements  corresponding  to  some  particular  solution.  Apply¬ 


ing  to  a  set  of  displacements  which  are  zero  at  the  boundary 
can  be  obtained. 


(2-12) 


The  system  [XQ]  can  be  considered  as  the' initial  nodal  forces' 
corresponding  to  the  surface  loading  on  the  single  shell.  These 
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are  the  nodal  forces  corresponding  to  zero  nodal  displacements 
under  the  action  of  the  applied  loading.  ^X1  J  can  also  be 
considered  as  a  new  particular  solution  with  zero  displace¬ 
ments  at  the  edges. 

The  stiffness  equations  including  consideration  of 
surface  loading  can  now  be  written  for  a  single  shell  element 
in  the  local  co-ordinate  system  as 

[X'l  =  [K]  [d']  +  [X'ol  •  (2-13) 

For  an  assembly  of  shells  the  total  structural  stiffness 
matrix  can  be  written  from  the  element  stiffness  matrices, 
transformed  to  global  co-ordinates  as  follows, 

[x]  =  [K]  [dl  +  pg  •  (2-14) 

In  EQUATIONS  (2-13)  and  (2-14)  [d1]  and  [dl  represent  the  net 
nodal  displacements. 

In  a  similar  manner  the  net  equivalent  nodal  forces  are 

given  by 

[Xe‘]  =  [K1]  [d1]  +  [x'0l  •  (2-15) 

From  the  equivalent  nodal  forces,  the  internal  forces  can  be 


obtai ned. 


( 


CHAPTER  III 


THEORY  OF  CLOSED  CIRCULAR  CYLINDERS  SUBJECTED  TO 
ARBITRARY  LOADING  ALONG  THE  EDGES 

3.1  INTRODUCTION 


In  this  chapter  the  exact  general  theory  proposed  by 
Flugge  is  outlined.  For  the  special  case  of  m  =  1 ,  the 
general  solution  given  by  Flugge  cannot  be  directly  applied 
because  of  the  repeated  roots  of  the  shell  characteristic 
equation.  The  solution  is  therefore  modified  for  m  =  1.  As 
an  exact  analysis  in  general  terms  is  tedious,  approximate 
methods  are  considered.  This  requires  a  prior  study  of  the 
different  shell  characteristic  equations  and  the  nature  of 
damping  of  the  solutions. 


3.2  GENERAL  THEORY 


The  exact  theory  has  been  proposed  by  Flugge  and  for 
a  complete  presentation  of  the  development  of  the  theory 
reference  may  be  made  to  pp  208  to  239, CHAPTER  5,  ref.  14. 


The  forces  and  displacements  (FIGURE  3.1)  are  periodic 
functions  of  <j>  and  hence  are  expressed  as  Fourier  series. 

Any  quantity  f,  either  a  force  or  a  displacement  is  expressed 
in  the  general  form  as 


f  =  c  [C  i  f  i  +  ^3^3 


+  C,f,  +  C,f,  +  Ccfc  +  Ccfc  +  C,f,  +  Cgfg] 


'4  4 


'5  5 


o  6 


7  7 


cos 

sin 


m<j)  (3-1 ) 
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The  constants 
four  at  each  edge, 
fol 1 ows* 

fi =  Vi  ■  V2 

f2  =  a2*  1  +  V2 

^3  —  ^3^2  ~  ^4^4 


to  Cg  are  evaluated 
The  quantities  f-j 


from  the  boundary  conditions 
to  fg  are  expressed  as 

^5  ^5^3  "*  a5$6 

f6  =  a6(p  5  ~  a5'6 

f7  ~~a7$ 7  _  a8$8 

fg  =  ag^7  -  ay$8  (3-2) 


The  characteristic  functions  $-j  to  $g  are  defined  as 


-y 

^  =  e  1  cosy-jy 

+x1y 

$5  =  e  cosy-jy 

-v  . 

+xxy 

$2  =  e  siny-jy 

t6  =  e  siny-jy 

-X2y 

+X2y 

$3  =  e  cosy2y 

$7  =  e  cosy2y 

-XgY. 

$4  =  e  siny2y 

+x  y 

$8  =  e  ^  siny2y 

(3-3) 

By  substituting  the  displacements  expressed  as  periodic 
functions  into  the  three  differential  equations  (EQUATIONS  13a-c, 
p  219,  ref.  14)  Flugge  obtaines  an  eighth  order  characteristic 
equation  (EQUATION  1-12).  The  roots  of  this  equation  occur  as 
pairs  of  complex  conjugate  quantities  expressed  in  terms  of 
x-j ,  y-j  and  x£>  ^2*  Values  °f  a]  ar|d  a2  have  been  Preser|ted  in 
a  tabular  form  (TABLE  1,  p.  230,  ref.  14). Values  of  a^  and  a^ 
can  be  obtained  in  turn  from  a-j  and  a2  as  shown  in  TABLE  3.1. 
Values  of  ag,a^,  a^  and  ag  are  obtained  using  these  two  tables 
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by  changing  the  subscripts  of  x,u>  from  1  to  2  and  those  of 
a, 3  from  1  and  2  to  3  and  4  respecti vely .  Knowing  the  shell 
parameters  and  harmonic  number  m  of  the  loading  the  roots  can 
be  obtained  from  the  characteristic  equation.  Values  of  a-j 
to  ag  can  be  determined  from  the  above  tables.  Using  these, 
values  of  f^  to  fg  can  be  found  at  any  point  along  the  length 
of  the  shell  by  means  of  EQUATIONS  (3-2)  and  (3-3).  Sub¬ 
stituting  these  into  EQUATION  (3-1),  the  value  of  any  stress 
resultant  or  displacement  f  at  any  point  along  the  length 
of  the  shell  can  be  expressed  in  terms  of  the  arbitrary 
constants  to  Cg. 

Since  in  the  derivation  of  EQUATION  (3-1)  the  stress 
resultants  and  other  displacements  are  expressed  in  terms  of 
the  displacement  w  it  is  pertinent  to  base  further  discussion 
on  w.  Using  TABLE  3.1  the  expression  for  w  can  be  written  as 

w  =  [e"xly[C'1cosy-1y  +  CgSin^y]  + 

e"x2y[C'3cosy2y  +  r4sinu2y]  +  e+xly[C5cosy-jy  + 

CgSin^y]  +  e+x2y[C7cosu2y  +  CgSiny2y]]cos  m<£  (3-4) 

There  are  two  cases  in  which  the  solutions  can  be 
simplified  within  the  framework  of  the  exact  theory.  These 
are  the  semi -infinite  cylinder  and  the  case  for  which  m  =  o. 

In  a  semi -infinite  cylinder  the  solutions  containing 
the  terms  e+xly  and  e+x2y  become  unbounded  as  y  increases.  For 


a  long  shell  if  the  applied  loads  are  in  self  equilibrium,  loads 
at  one  edge  will  not  produce  stresses  and  displacements  at  the 
other  edge.  We  can  therefore  disregard  the  four  solutions 
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correspond]' ng  to  to  Cg.  The  expression  for  w  can  be  written  as 


e  x2y[C3cosy-|y  +  C^siny-|y]  cos  mf 


(3-5) 


The  case  m  =  o  gives  rise  to  axi symmetric  bending 
which  has  been  fully  covered  in  literature.  This  case  is 
governed  by  a  fourth  order  differential  equation  (EQUATION  73a, 
p.  273,  ref.  14)  giving  rise  to  a  solution  in  terms  of  four 
arbitrary  constants  of  integration.  The  expression  for  w 
can  be  wri tten  as  > 


w  =  e_M^[C-j  cosyy  +  C^sinyy]  + 
e+yy[C3cos“y  +  C4sinyy] 


(3-6) 


As  before,  ignoring  the  interaction  of  the  far  edge  for  a 


long  shell,  in  this  expression  the  second  term  involving  e+yy 


can  be  disregarded,  thus  reducing  the  number  of  constants  to  two. 

The  case  m  =  1  requires  special  treatment  in  that  the 
values  of  x2>  ^2  are  zero*  As  the  f°ur  roots  expressed  in  terms 
of  x2>  u2  are  repeated  and  zero  the  part  of  the  general  expression 
for  w  given  by  EQUATION  (3-4)  is  not  valid  and  must  be  re¬ 
placed  by  a  simple  polynomial  solution.  The  complete  solution 
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then  takes  the  form, 

=  [e  x1y[C1cosy1y  +  C^sin^y]  +  e+xl y[C’5cosy + 

cos  m 6 


w 


CgSin^y]  +[C3  +  C4y  +  C?y2  +  Cgy3] 


TABLE  3.2  gives  the  displacements  and  stress  resultants 
corresponding  to  the  polynomial  part  of  the  solution.  Using 
this  table,  the  membrane  stresses  can  be  written  as, 


Nx  =  D  (1-u2)[2C1  +  6yCg] 

V  v =  6“  (1-u2)  s 


(3-7) 


(3-8) 


N  ,  =  o 

<P 

This  table  gives  also  moments  and  transverse  forces.  These 
quantities  are  of  the  order  k  and  therefore  the  resulting 
stresses  are  insignificantly  small  in  comparison  with  those 
produced  by  membrane  forces.  Thus  they  can  be  ignored. 
Expressions  identical  to  EQUATIONS  (3-8)  can  also  be  obtained 
from  the  membrane  theory  for  m  =  1.  Therefore  for  a  thin 
shell  the  polynomial  part  of  the  general  solution  can  be 
exactly  replaced  by  the  membrane  solution. 


It  is  seen  that  for  the  particular  cases  of  m  =  o 
and  a  long  shell  the  number  of  constants  to  be  determined 
from  the  boundary  conditions  are  reduced.  In  a  general  case 
one  has  to  evaluate  all  the  eight  arbitrary  constants.  The 
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evaluation  of  these  constants  in  general  terms  using  an 
exact  theory  as  described  above  is  very  tedious  and  time 
consuming.  However  values  of  stress  resultants  can  be 
tabulated  for  particular  cases  of  T/r,  h/r,  m  and  u.  A 
table  similar  in  form  to  the  one  given  in  A.S.C.E.  Manual  31 
for  open  cylindrical  shells  for  loads  applied  at  the  straight 
edges  can  be  obtained.  Such  tables  would  have  to  be  very 
exhaustive  to  consider  sufficient  values  of  these  parameters 
to  include  all  practical  cases.  As  an  example  TABLE  3.3  has 
been  presented  giving  values  of  the  stress  resultants  form  =  1 
and  2,  for  particular  values  of  the  shell  parameters. 

A  solution  in  closed  form  in  general  terms  may  be 
possible  if  one  uses  approximate  methods  to  the  general  shell 
problem.  This  requires  a  study  of  the  different  shell 
characteristic  equations  and  the  nature  of  the  solutions  and 
this  is  considered  in  the  next  section. 

3.3  STUDY  OF  CHARACTERISTIC  EQUATIONS 

Flugge's  exact  characteristic  equation  and  Holand's 
modified  equation  are  used  for  an  accurate  analysis  of  the  shell 
problem.  Donnell's  characteristic  equation  is  widely  used 
because  of  its  relative  simplifity.  Two  independent  fourth 
order  equations  in  place  of  a  single  eighth  order  equation  have 
been  proposed  by  Fl'ugge  for  low  values  of  m.  (refer  CHAPTER  I). 
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The  explicit  solutions  for  the  roots  of  Hoi  and,  Donnell  and 
the  split  fourth  order  equations  of  Flugge  have  been  given  in 
APPENDIX  A. 

Flugge 's  eighth  order  characteristic  equation  is 
considered  exact.  The  solutions  of  the  other  equations  are 
therefore  compared  with  those  of  Flugge.  A  value  of  0.167 
has  been  adopted  for  Poisson's  ratio  u,  this  being  the  usual 
value  for  concrete.  It  has  been  found  that  variation  in  u 
does  not  affect  the  roots  appreciably.  Therefore  the 
qualitative  nature  of  any  further  discussion  in  this  thesis 
is  not  affected  by  variation  in  o. 

FIGURES  3.2  and  3.3  show  the  variation  of  the  quantities 
X-j  ,  y i  and  X2>  u2  for  various  values  of  r/h  computed  using 
the  exact  theory.  It  is  seen  from  the  figures  that  values  of 
x-j  and  x2  approach  y-j  and  respectively  as  r/h  increases.  This 
property  is  more  pronounced  for  lower  values  of  m. 

TABLE  (3-4)  show  a  comparison  of  the  roots  obtained 
from  the  different  characteri stic  equations.  Holand's  roots 
are  found  to  agree  with  those  of  Flugge  for  all  values  of  m. 

The  agreement  is  quite  good  as  k  becomes  very  small.  For  values  of 
m  considered  in  this  table  xy  y-j  obtained  by  Donnell's  equation 
have  a  better  agreement  with  those  of  Flugge  than  the  ones 
given  by  Holand.  For  m  =  1,  Donnell's  equation  does  not  give 
the  zero  values  of  x2 »  ^  which  are  essentia1  for  the  problem. 
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For  lower  values  of  m,  values  of  x£>  obtained  by  Donnell's 
equation  do  not  agree  with  the  exact  values.  However  the 
agreement  improves  as  m  increases,  provided  m  is  not  very  high. 
For  lower  values  of  m  the  roots  obtained  by  the  splic  equations 
of  Flligge  agree  well  with  the  exact  values.  The  agreement  be¬ 
comes  better  as  k  becomes  smaller.  For  m  =  1,  values  of  x£> 
are  zero  indicating  a  perfect  agreement  with  the  exact  values. 

3.4  STUDY  OF  THE  NATURE  OF  DAMPING  OF  THE  SOLUTIONS 


The  nature  of  approximations  to  the  shell  problem  depends 
on  the  nature  of  damping  of  the  solutions.  Therefore  this  aspect 
is  investigated  in  detail  in  this  section. 


Making  the  substitutions  y=Q-y  in  the  last  two  terms 
of  EQUATION  (3-4)  we  get  the  expression  for  w  as 


e"xly[C1cosu1y  +  C^ini^y]  +  e“x2-y  [C3cosu2y  +  C4siny2y]  + 


w  = 

rx-,7— ' 


“XoY rr 


—  —  i  _ 


e"xly[C5'cosy-1y  +  Cjsiny^]  +  e~x2y[C7'cosy2y  +  C8siny2y] 

where  the  constants  CQ-  to  Cg  are  modified  to  Cg  to  Cg  to  take 
into  account  this  transformation. 


cosine})  (3-9) 


EQUATION  (3-9)  now  consists  of  four  different  components  each 
representing  damped  oscillations  like  the  one  shown  in  FIGURE  3.4. 

As  cosily  is  periodic,  the  period  of  oscillation  lT  is  given  by 
the  expression, 
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L 


2tt 

y 


where  i  =  1,2 


As  y  is  nondimensional ,  the  absolute  value  of  this  period,  L. 


can  be  written  as, 


L. 

i 


2irr 


U 


l 


L  represents  the  length  of  the  shell  over  which  any  edge  distur- 

i 

bance  again  reaches  a  maximum  value.  The  ratio  of  these  maximum 
values  represents  the  rate  of  damping  which»referring  to  FIGURE  3.4 
is  given  by  the  expression, 

A 


2ttX, 


e  xi^ 


exi  L . 

i  _  e 


y  i 


A 


n+l 


e_xi  (y+  L  •) 


(3-11) 


For  low  values  of  m  it  is  seen  from  FIGURES  3.2  and  3.4 
that  y  .  is  approximately  equal  to  Xj  .  Therefore  EQUATION  (3-11)  can 
be  written  as. 


A. 


2tt 


^n+1 


=  e 


=  535.5 


Thus  for  low  values  of  m,  the  rate  of  damping  is  independent  of  the 
dimensions  of  the  shell. 


For  high  values  of  m  asy^is  always  less  than  Xi  the  rate  of 

2tt 

damping  is  greater  than  e  .  This  means  that  for  all  values  of  m 


the  rate  of  damping  is  very  quick. 


The  purpose  of  the  above  study  is  to  determine  when  a  distur- 
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bance  applied  at  the  near  edge  may  be  neglected  at  the  far  edge. 

EQUATION  (3-10)  indicates  that  the  period  is  a  function  of  the 

shell  dimensions  and  the  harmonic  number  m.For  example  TABLE  3.5  gives 

values  of  the  periods  U|  and  {A,  f°r  a  moderately  thick  and  a 

thin  shell  for  two  different  harmonics. Consider  first  the 

low  harmonic  m  =  2.  For  the  thick  shell  for  which  h/r  =  0.03, 

taking  l/r  =  1,  the  x-j-solution  is  quickly  damped,  but  the 

solution  is  undamped.  However  when  the  length  is  very  large 

say  Z/r  =  25,  the  X2's0^ut’1'0n  is  a^so  damped.  For  thinner  shells 

the  situation  is  slightly  different.  For  the  thinner  shell  for 

which  h_  =  0.0014  ,  for  all  the  lengths  the  x-|- solutions  are 

r 

quickly  damped,  but  the  X2  “  solutions  are  practically  un¬ 
damped.  For  the  higher  harmonic  m  =  10  the  length  over  which 
damping  occurs  is  slightly  increased  for  the  Xj  solution,  but  this 
length  is  considerably  reduced  for  the  X2  solution. 

From  a  study  of  the  roots  of  the  characteristic 
equations  and  the  nature  of  damping  of  the  solution  the 
following  concl usions  can  be  drawn. 

1.  Holand'S  roots  can  be  expected  to  give  results 
which  are  pretty  close  to  the  exact  values  obtained  from  Flugge. 

2.  As  the  x-r  solution  is  quickly  damped  for  all 
values  of  &/r,  h/r  and  m,  slight  variation  in  the  values  of 
X-j ,  p-j  cannot  affect  the  final  solutions  appreciably. 

Incidentally  this  means  that  for  all  values  of  m,  except  very  high 
values  Donnell's  equation  can  be  used  to  obtain  the  solution 


V 


‘  t 


corresponding  to  Xj  J  ^i* 

3.  As  the  x£" sol ution  is  slowly  damped  for  low  values 

of  m,  slight  variation  in  the  values  of  X2>  ^2  cdn  a^^eci: 
the  final  solutions  appreciably.  For  low  values  of  m  the 
values  of  are  eva^ua'te<^  with  a  higher  degree 

of  accuracy.  This  means  that  for  low  values  of  in,  Donnell's 
equation  cannot  be  used  to  obtain  the  solution  corresponding 
to  X£  *  r  2 • 

4.  As  for  low  values  of  m,  values  of  Xy  r-j  andx£, 

computed  by  the  split  equation  of  Flugge  are  close  to  the 
exact  ones,  these  equations  can  be  expected  to  give  comparabl 


results . 
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TABLE  3.2 

CHARACTERISTIC  COEFFICIENTS  FOR  THE 
POLYNOMIAL  PART  OF  THE  SOLUTION 


f 

c 

f3 

f4  f7 

f8 

¥ 

1 

1 

y 

2 

y 

3 

y 

¥' 

1 

0 

1 

2y 

o  2 

3y 

u 

1 

0 

-1 

-2y 

-6(2+u)-3y2 

V 

1 

-1 

-y 

2u-y2 

6uy-y3 

D 

0 

0 

0 

0 

"x 

D 

0 

0 

2[u2-(Hk)] 

[u2-(l+k)  ]  6y 

V 

D 

0 

0 

0 

6(l-u2)+3k(2+u) 

% 

D 

0 

0 

0 

6(l-u2)+3k(u-u2) 

Tx 

D 

0 

0 

0 

6(1 -o2)+9k(u-u2) 

Dk 

0 

0 

2u 

6uy 

M 

X 

Dk 

0 

0 

2(2-u2) 

6(2-u2)y 

M~ 

(J)X 

Dk 

0 

0 

0 

6(1 -o) 

Dk 

0 

0 

0 

-6u(l -u) 

% 

Dk 

0 

0 

-2u 

-6oy 

«x 

Dk 

0 

0 

0 

6[3-u-u2] 

0 

Ok 

0 

0 

0 

6(3-2u) 

f 

=  c 

^3^3  ^4^4  *  ^7 

f7  + 

38 


co 

CO 

LiJ 

_J 

CQ 

f— 


00 
UJ 
ZD 
— I 
< 
> 


o 


oo 

h- 


< 

h- 


oo 

UJ 

c£ 

oo 

oo 

LU 

CrL 

I— 

oo 


TABLE  3.4 


VALUES  OF  ROOTS  OBTAINED  BY  FLUGGE,  HOLAND 
DONNELL  AND  APPROXIMATE  SOLUTIONS 
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k  m 

X1 

P1 

X2 

VZ 

0.8333x1 0"3  1 

Flligge 

4.2419 

4.0200 

0.0 

0.0 

Hoi  and 

4.1325 

4.1325 

0.0 

0.0 

Donel 1 

4.2568 

4.0152 

0.1243 

0.1173 

Approx* 

4.1325 

4.1325 

0.0 

0.0 

5 

Fliigge 

7.1238 

2.8280 

2.9407 

1.2513 

Hoi  and 
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2.9164 

2.9555 

1.2161 
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7.1366 
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3.0042 

1.2242 
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5.8754 
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2.0845 
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*  Approx. 

Rased  on  split  equations  of  Fliigge. 
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TABLE  3.5 
PERIODS  OF  DAMPING 
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2 
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126.00 

0.0014 

10 
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FIGURE  3.2  VAR f AT J OKI  OF  BY  EXACT  THEORY 


43 


FIGURE  3.3  VARIATION  OF  Xz  }/Uz  BY  EXACT  THEORY 
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EfGURE  3-4  MATURE  OR  PA  MP/A/G> 


CHAPTER  IV 


INFLUENCE  COEFFICIENTS 

4.1  INTRODUCTION: 

In  this  chapter,  the  method  of  finding  the  stiffness  in- 
fluence  coefficients  for  the  closed  cylindrical  shell  using 
the  exact  theory  is  described.  A  specific  numerical  example 
has  been  given  to  illustrate  the  method.  The  resulting  direct 
stiffness  matrix  contains  64  elements  and  can  be  made  truely 
symmetric  by  adopting  the  stiffness  matrix  sign  convention. 

However  by  adopting  a  slightly  different  sign  convention  this 
matrix  can  be  reduced  such  that  only  20  elements  are  different. 

Establishing  the  influence  coefficients  in  general  terms 
is  tedious.  Therefore  approximate  methods  are  attempted  and 
these  are  discussed  in  detail.  It  is  seen  that  three  types  of 
approximations  are  possible  depending  on  the  shell  geometry  and 
the  harmonic  number  of  the  applied  edge  displacement.  A  criterion 
for  the  limitation  of  replacing  the  basic  solution  by  the  membrane 
solution  has  been  established. 

4.2  METHOD  OF  FINDING  INFLUENCE  COEFFICIENTS  BY  THE  EXACT  THEORY 

For  any  particular  harmonic  m  the  four  displacements  w',  w, 
u  and  v  at  each  of  the  two  edges  of  a  closed  cylinder  can  be 
expressed  in  terms  of  the  eight  constants  C-j  to  Cg  as  explained 
in  CHAPTER  III.  Unit  value  of  each  displacement  is  applied  individually 
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at  each  edge  while  the  remainder  of  the  displacements  are  set  as 
zero.  The  resulting  equations  obtained  from  EQUATION  (3-1)  can 
be  written  in  matrix  notation, with  [i ]  as  the  unit  matrix, as 

[C-,]  [C]  =  [I]  (4-1) 

The  coefficient  matrix  [C^]  of  size  8x8  found  from  the  elements  f-j 
to  fg  is  independent  of  the  boundary  conditions.  Each  column 
vector  of  the  matrix  [c]  consists  of  the  unknown  constants  C-| 
to  Cg  corresponding  to  each  column  vector  of  [I].  The  elements 
in  each  of  this  column  vector  represents  a  set  of  displacement 
boundary  conditions  for  each  case  of  the  applied  unit  edge  dis¬ 
placement.  In  the  analysis  that  follows,  the  displacements  have 
been  considered  in  the  order  w',  w,  u  and  v  at  each  edge. 

Solving  EQUATION  (4-1)  we  get 

[C]  =  [qrki]  =  [Cjf1  (4-2) 

The  direct  stiffness  matrix  can  now  be  determined.  The 

stress  resultants  M  ,  S  ,  N  and  T  considered  in  this  order  to 

xxx  x 

correspond  with  the  order  of  displacements  can  be  expressed  in 
terms  of  the  eight  constants  C-j  to  Cg  for  each  case  of  the 
applied  edge  displacements.  The  stiffness  matrix  can  be  written 
as , 

-1 


[K] 


[c2]  [eg 


(4-3) 
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From  EQUATION  (3-1)  in  a  manner  similar  to  that  of  [C-j]  for 
displacements,  the  coefficient  matrix  [C^]  for  the  stress 
resultants  is  established,  and  as  such  is  independent  of  the 
boundary  conditions. 

The  crux  of  the  problem  lies  in  evaluating  the  constants 
given  by  each  column  vector  of  [C].  This  requires  determination 
of  the  inverse  of  [C-j  ] .  EQUATION  (4-4)  gives  the  matrix  [C-j] 
in  general  terms  of  shell  geometry  and  loading.  The  inverse 
can  be  readily  obtained  if  specific  numerical  values  are  used 
for  these  terms.  In  the  course  of  the  study  the  inverse  of 
[C  -j  ]  was  in  fact  obtained  in  general  terms,  but  the  complexity 
of  the  resulting  expressions  was  such  that  they  were  not  suit¬ 
able  for  practical  applications. 

An  example  illustrating  the  application  of  the  exact 
theory  as  outlined  above  follows.  EQUATION  (4-5)  gives  the 
stiffness  matrix  [K-j]  for  a  particular  case  of  Q=2,  P=0.03, 
u  -0.1 67  and  m=2.  The  matrix  was  evaluated  using  a  digital  computer. 
The  sign  convention  adopted  is  as  shown  in  FIGURE  4.1  which 
is  based  on  the  sign  convention  used  in  the  general  shell  theory 
(FIGURE  3.1).  Since  the  parameters  Q  and  P  are  non  dimensional, 
to  obtain  the  absolute  values  of  the  elements  in  the  stiffness 
matrix  for  absolute  values  of  the  edge  displacements  all  the 
elements  in  the  1st  and  the  5th  rows  and  all  the  elements  in 
the  1st  and  the  5th  columns  are  to  be  multiplied  by  r. 
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It  is  noted  that  the  stiffness  matrix  [K^]  in  EQUATION  (4-5) 
is  not  symmetrical  with  regard  to  signs.  This  results  because 
of  the  sign  convention  used  in  the  general  shell  theory.  A 
truly  symmetrical  stiffness  matrix  can  be  obtained  by  adopting 
the  usual  sign  convention  for  the  stiffness  method.  This  sign 
convention  is  given  in  FIGURE  4.2  and  the  stiffness  matrix 
[K^]  using  this  sign  convention  is  given  in  EQUATION  (4-6).  In 
accordance  with  this  classical  stiffness  matrix  formulation  [ ] 
can  be  written  as  follows, 


(4-7) 


B 


C 


In  the  matrix  [ ]  the  numerical  values  of  the  elements 
of  [C]  and  [B^]  are  identical  to  those  of  [A]  and  [B]  respectively. 
It  would  be  convenient  to  introduce  a  sign  convention  which  would 
also  make  the  algebraic  signs  identical  so  that  the  matrix  can 
be  written  as  in  EQUATION  (4-8)  * 


A  B 

B  A 


(4-8) 


The  sign  convention  for  this  matrix  is  given  in  FIGURE  4.3  and  the 
matrix  [Kg]  for  the  above  problem  written  using  this  sign  con¬ 
vention  is  given  in  EQUATION  (4-9).  It  is  seen  that  the  matrix 
[Kg]  is  not  only  symmetrical  but  also  that  it  can  be  formed  by 
two  symmetrical  submatrices  [A]  and  [B].  This  means  that  out 
of  the  64  elements  of  the  matrix  [Kg]  only  20  elements  are  different. 
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As  explained  earlier  it  is  very  difficult  to  obtain  the 
stiffness  influence  coefficients  in  general  terms  because  of  the 
difficulty  in  determining  the  inverse  of  [C-j].  Since  one  of 
objectives  of  this  study  is  to  obtain  general  expressions  for 
stiffness  influence  coefficients  approximate  solutions  are 
consi dered. 

4.3  APPROXIMATE  SOLUTIONS 

In  accordance  with  Goldenveiser  and  Flligge  as  discussed 
in  the  literature  review  in  CHAPTER  1,  the  approximate  solutions 
to  the  problem  of  higher  harmonic  loading  on  the  edges  of  closed 
cylinders  are  as  outlined  below: 

1.  The  split  type  solution  for  low  values  of  m,  where 
the  tangential  and  non  tangential  boundary  conditions  are  con- 
si  dered  separately. 

2.  The  general  solution  for  high  values  of  m,  where  both 
the  tangential  and  non  tangential  boundary  conditions  at  one 

edge  are  considered  together  neglecting  the  effect  of  the  far  edge. 

3.  The  solution  for  very  high  values  of  m,  where  the 
problem  is  split  into  two  states  of  stress  corresponding  to  general¬ 
ised  plane  stress  and  plate  bending. 

There  are  thus  3  types  of  approximations  possible  depending 
on  the  loading  harmonic,  and  stiffness  influence  coefficients  based 
on  these  approximations  are  presented  in  this  chapter. 
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4 . 4  BASED  ON  SPLIT  TYPE  SOLUTION  FOR  LOW  VALUES  OF  m 

The  general  solution  is  split  into  two  parts,  the  edge 
effects  and  the  basic  state.  The  edge  effects  part  is  governed 
by  the  characteristic  equation, 


X4  -  2 (2m2  -  u)x2  +  [1  -  u2  +  6m2 (m2  -  1)]  =  o 

k 


(4-10) 


and  the  basic  state  part  by  the  equation, 


[1  -  u2  +  6m2 (m2  -  1)]X4  -  2m2 [2m4  -  (4  -  u)m2 


k 

+  (2-u)]  X  +  m4(m2  -  l)2  =  o 


(4-11) 


These  equations  have  been  presented  previously  in  CHAPTER  1  as 
EQUATIONS  ( 1  - 1 3a , b ) .  The  equations  are  considered  separately  in 
the  following  sections. 

4.4.1  EDGE  EFFECTS  SOLUTION 

The  evaluation  of  the  roots  for  EQUATION  (4-10)  is  given 
in  SECTION  A. 3  of  APPENDIX  A.  Making  the  assumption  that, 


1  -  u2  +  6m2 (m2  -  1 )  >>  (2m2  -  u) 
k 


(4-12) 


we  obtain  > 


-4  -7 

Since  k  is  very  small  and  generally  of  the  order  10  to  10 
and  as  m  is  also  low  this  assumption  is  justified. 


(4-13) 


55 


The  equality  given  by  EQUATION  (4-13)  leads  to  further 
simplifications.  Making  also  the  assumption  that 


analytical  expressions  for  a-j ,  c^,  3-j ,  3^  have  been  derived  in 
SECTION  A. 5  of  APPENDIX  A.  These  terms  relate  the  three  displacements 
u,  v  and  w  and  a  complete  explanation  of  these  terms  is  presented 
on  p.  228,  Chapter  5,  ref.  14.  For  the  latter  assumption  to  be 
valid  x-j  must  be  much  greater  than  m,  a  condition  that  is 
satisfied  for  small  values  of  k. 

Values  of  a-j ,  a^,  3-|  ,  3 ^  bave  been  evaluated  for  various 
values  of  m  and  k  and  presented  in  TABLE  4.1.  The  values 
labelled  Flligge  and  Donnell  in  this  table  are  obtained  using 
exact  solutions  to  their  respective  equations  without  any  approxi 
T-mations  whereas  the  simplified  solutions  are  based  on  the  split 
equations  of  Flugge  with  the  above  two  assumptions. 

From  the  table  it  is  seen  that  the  values  of  a-,,  3-j 

and  32  found  from  the  simplified  solutions  agree  well  with  the 
exact  values  provided  m  and  k  are  very  small. 

TABLE  4.1  also  illustrates  that  for  all  values  of  m  and  k 
values  of  the  parameters  a.  and  3.  obtained  from  x]  based  on 
Donnell's  theory  agree  well  with  the  exact  values  of  Flugge,  in¬ 
dicating  that  the  Donnell's  theory  can  be  used  for  this  part  of 
of  the  solution  as  concluded  in  SECTION  3.4. 
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For  the  split  form  of  the  EQUATION  (4-10),  the  non 
tangential  boundary  conditions  corresponding  to  w 1  and  w  are 
satisfactory  for  establishing  the  influence  coefficients  for 
the  edge  effects  solution.  Further  if  x-j  =  y-j  ,  the  solution 
takes  the  form* 

w  =  C-j(e"xly  cos  x-jY) +  C2(e  ~xly  sin  x-jY) 

+  C5(e  +xly  cos  X]y)+  Cg(e  +xly  sin  x-jY)  *  (4-15) 

As  shown  by  the  values  given  in  TABLE  3.4,x^  is  generally  large. 
Therefore  the  last  two  terms  become  unbounded  as  y  increases 
while  the  first  two  terms  are  highly  damped.  Because  of  the 
nature  of  the  x-|~  solutions  the  last  two  terms  can  be  neglected 
ignoring  the  effect  of  the  far  edge.  The  solution  can  now  be 
wri tten  as , 

w  =  C'-jC e“xly  cos  x-|Y)+  C‘2(e_xly  sin  x-|Y)  '  (4-16) 

Designating  e"xly  cos  x-j Y  and  e~xr  sin  x-jY  as  g  and  g^res- 

_  I 

pectively,we  can  write  w  and  w  as  , 
w  -  ^  ^2^2  9 

w‘  =  C-j  ( —  x 1 9 1  -  x 1 9 2 )  +  ^2^"xlg2  +  xlgl  )  ~ 
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For  the  boundary  conditions  w1  =  1  and  w  =  o  at 
the  top  edge  the  following  equations  are  obtained  for  the 
constants  C-j  and  C^, 


“X1 

X1 

cl~ 

1 

) 

1 

0 

>C2- 

0 

from  which, 

C-j  =  o  and  =  1_  * 

X1 

Using  a  relation  similar  to  EQUATION  (3-1)  any  quantity  f, 
either  a  stress  resultant  or  a  displacement,can  be  expressed  as, 

f  =  c  [( a i C*-j  +  a2C2)g1  +  (a^  -  a^jgj  . 

Substituting  values  of  C-j  and  we  get, 

f  =  C  [^2  g,  +  !i  g2l  •  (4-17) 

X1  X1 

Values  of  a^  and  a2  for  all  the  displacements  and  stress 
resultants  are  given  in  TABLE  4.3.  This  table  is  obtained  in 
the  following  manner.  Using  the  simplified  expressions  of  Flugge 
for  the  stress  resultants  (see  EQUATIONS  (12a-f),  p.  217,  ref.  14), 
the  table  of  characteristic  coefficients  (TABLE  1,  p.  230,  ref.  14) 
has  been  reduced  to  a  much  simpler  form  as  given  by  TABLE  4.2. 

Making  the  substitution  =  x-j  in  this  table  we  obtain  TABLE  4.3. 
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TABLE  4.4  which  gives  the  characteri Stic  coefficients  for  the 
boundary  conditions  w1  =1  and  w  =  o  at  the  top  edge  is 
obtained  by  substituting  the  values  of  a-j  and  a2  from  TABLE  4.3 

in  EQUATION  (4-17).  Using  TABLE  4.4  explicit  expressions  for 

«««• 

any  displacement  or  stress  resultant  f  for  these  boundary 
conditions  can  be  written  as, 

f  =  c^S-j  g-j  +  S2g2)  (4-18) 

For  the  boundary  conditions  w'  =  o  and  w  =  1  at  the 
top  edge  the  following  equations  are  obtained  for  the  constants, 


~"xl 

xl' 

"  f  1 
L1 

0 

_  1 

0 

_c2. 

1 

from  which 

C-j  =  1  and  C2  =  1 . 

As  before  any  quantity  f  can  be  expressed  as 
f  ~  c[  (a i C -j  +  a 2C 2 ) 3 1  ^  (a] ^2  ~  a2^l  ^2^ 

Substituting  the  values  of  C-j  and  C2  and  a^  and  a2  from  TABLE  4.3^ 

f  =  c  (  s3g1  +  s^g2  )  (4-19) 

In  a  similar  manner  TABLE  4.5  giving  the  characteristic 
coefficients  for  the  boundary  conditions  w'  =  o  and  w.  =  1 
at  the  top  edge  is  obtained.  Using  this  table  explicit  expressions 
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for  stress  resultants  and  displacements  for  these  boundary 
conditions  can  be  written. 

From  the  expressions  given  in  TABLES  4.4  and  4.5,  the 

displacements  w' ,  w,  u  and  v  and  the  stress  resultants  M  ,  S  , 

x  x 

N  and  T  at  the  top  edge  correspondi ng  to  y  =  o  can  be  obtained. 

X  X 

Both  these  displacements  and  stress  resultants  are  quickly 
damped  because  of  the  nature  of  the  x-p solution.  The  values 
at  y  =  o  are  obtained  by  making  g-j  =  1  and  =  o  in  the 
expressions  given  in  TABLES  4.4  and  4.5. 


It  should  be  pointed  out  that  as  the  other  boundary 
conditions  u  =  o  and  v  =  o  have  not  been  considered  these  quantities 
will  appear  as  non  zero  displacements  at  the  top  edge.  The 
effect  of  these  displacements  is  discussed  later  in  SECTION  4.4.3. 


The  stiffness  matrix  correspondi ng  to  the  edge  effect 


solution  can  be  written  from  the  values  at  y  =  o.  As  m  «  x-| 

2 

in  the  expressions  for  M  and  S  the  terms  involving  m  can  be 

X  X 

2 

neglected  compared  to  x-j  •  A  symmetric  matrix  [SM]]  is  then 


as  shown 

bel ow. 

w  '  =  1 

w  =  1 

[SMi  ]  =  K 

"lxl 

-v 

_-2Xl2 

8x,3  . 

M 


x 


x 


(4-20) 
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where  the  constant  multi  pier  K  in  this  expression  is  the  shell 
flexural  stiffness. 

4.4.2  BASIC  SOLUTION 

This  solution  is  obtained  from  the  characteristic 
EQUATION  (4-11).  The  X2  solution  of  this  equation  is  very  slowly 
damping.  It  is  found  from  investigations  by  the  author  that 
the  final  stress  resultants  are  highly  sensitive  to  the  values 

of  ^2’  a3’  a4’  ^3  anc*  ^4  anc*  even  sma^  variations  in  these 
values  will  cause  the  matrix  for  the  evaluation  of  the  arbitrary 
constants  from  the  boundary  conditions  to  be  highly  ill  conditioned. 
This  precludes  introducing  approximations  to  simplify  the  ex¬ 
pressions  to  determine  these  quantities.  Their  determination  by 
the  exact  theory  is  tedious.  However  the  values  of  X2»  ^  are 
very  small  and  almost  zero.  If  x2»  ^2  are  considered  as  zero 
the  basic  solution  given  by  EQUATION  (4-11)  reduces  to  the 
following  form 

w  =  (D  +  Cy  +  By^  +  Ay^)  cos  mf  (4-21) 

This  solution  is  identical  to  the  one  obtained  by  the 
membrane  theory  (SECTION  A. 8  of  APPENDIX  A).  Hence  the  basic  solution 
can  be  replaced  by  the  membrane  theory  as  a  good  approximation. 

The  limitation  of  the  assumption  that  the  basic  solution 
can  be  replaced  by  the  membrane  solution,  can  be  established  in 
the  following  manner.  The  term  e 


in  the  basic  solution 
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corresponding  to  EQUATION  (4-11),  can  be  written  as, 

2  2  3  3  4  4 

Xy  AY  X  y  X  y  X  y  yXy 

l!  2!  3!  4! 

where  y  is  such  that  o<y<l .  Taking  the  polynomial  solution 
given  by  EQUATION  (4-21)  implies  the  assumption 

X4  y4  ^  <<  1 

As  XX  y  is  very  small  eY^  y  -*  1 . 

To  obtain  the  largest  error  set  y  =  Q. 


4! 

It  can  be  seen  from  EQUATION  (4-11)  that  X  -  h^k^m2,  where  the 
modulus  of  h ^  is  of  the  order  unity.  Substituting  this  value  in 
the  above. 


i/2  \2  4  n4 

k(m  -  1 )  m  Q 

41 


<<  1  or 


n2n4,  2  ^2  4 

P  Q  (m  -  1 )  m 

1 2(41/ 


«  1 


(4-22) 


Let  us  designate  the  quantity  on  the  left  hand  side  of  this 
inequality  as  p.  p  sets  a  limit  on  replacing  the  basic  solution 
by  the  membrane  solution.  For  m  =  1  ,  it  is  seen  that  p  =  o 
indicating  that  the  basic  solution  is  exactly  replaced  by  the 


membrane  solution. 
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1 
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0 

6 

m2Q3-2(l+u)Q 

6 

mQ2 

mQ 

1 

D 

k  - 

L°  - 

(4-23) 
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Influence  coefficients  are  now  determined  using  the 
membrane  theory  for  edge  loading  as  given  in  SECTION  A. 8  of 
APPENDIX  A.  The  tangential  boundary  conditions  u  and  v  are 
considered  for  establishing  the  influence  coefficients. 

For  the  boundary  conditions  u  =  1  and  v  =  o  at  the 
top  edge  and  the  correspond!' ng  displacements  at  the  bottom 
edge  set  as  zero,  the  EQUATIONS  (4-23)  are  obtained  for  the 
arbitrary  constants  A,B,C  and  D.  Values  of  these  constants 
given  in  TABLE  4.6  are  obtained  from  the  solution  of  these 
equations . 

Any  displacement  or  stress  resultant  f  can  be  expressed  as, 

f  =  c(F4  +  F3y  +  F2y2  +  F^y3 )  •  (4-24a) 

The  expressions  for  F|  to  given  in  TABLE  4.7  can  be  obtained  by 
substituting  the  values  of  the  constants  A,B,C  and  D  into 
the  approximate  stress  displacement  relations  given  by  EQUATIONS 
1 2a-f ,  p.  217,  ref.  14- 

For  the  other  pair  of  boundary  conditions  u  =  o  and  v  =  1 
at  the  top  edge  and  the  correspond!* ng  displacements  at  the  bottom 
edge  set  as  zero,  EQUATION  (4-25)  is  obtained  for  the  arbitrary 
constants.  Values  of  the  constants  given  in  TABLE  4.6  were 
obtained  from  the  solution  of  this  equation. 
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2n 

n  Q 
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"I 

Q 


n 

1 


-n_ 

2 


b_ 

Q 


1  N 


Q 


t 

Xcf> 


-n  1  -n 

2  Q  2 


Nx*  (4-26) 


2  2 

Note:  b  =  m  Q  /3  +  2(l+u) 
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As  before  the  corresponding  displacements  and  stress 
resultants  can  be  expressed  as 

f  =  c  (  Fg  +  F?y  +  F6y2  +  F^3  )  (4-24b) 

where  the  values  of  F-j  to  are  given  in  TABLE  4.8. 

In  both  cases,  the  membrane  stress  resultants  are  obtained 
by  considering  only  the  constants  A  and  B.  The  constants  C  and  D 
depend  only  on  the  inextensional  deformation  and  give  rise  to 
the  stress  couples.  For  m  =  1  ,  the  inextensional  deformation 
is  replaced  by  rigid  body  deformations  and  therefore  C  and  D  can 
be  ignored. 

From  the  expressions  given  in  TABLES  4.7  and  4.8  the 
stress  resultants  M  ,  S  ,  N  and  T  at  the  top  and  bottom  edges 
corresponding  to  y  =  o  and  y  =  Q  respectively  can  be  obtained. 

The  stiffness  matrix  corresponding  to  the  basic  part  can  now 
be  written  using  these  values.  A  symmetric  matrix  [ SM? ]  is 
obtained  as  shown  in  EQUATION  (4-26). 

4.4.3  COMPLETE  STIFFNESS  MATRIX 

To  obtain  a  complete  stiffness  matrix  for  a  shell  element 
it  is  necessary  to  assign  in  turn  to  each  boundary  displacement 
a  value  of  unity,  while  all  other  displacements  are  kept  zero. 
Therefore  the  elements  obtained  for  the  submatrices  [SMl]  and 
[SM2]  in  the  previous  sections  cannot  be  elements  in  the  complete 
stiffness  matrix  since  the  interaction  of  the  two  systems  of 
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boundary  condi tions.  corresponding  to  .the  edge. effects  and  , the 
basic  state  has  not  been  considered  so  far.  This  results  in  edge 
displacements  u  and  v  for  the  edge  effects  state  and  edge  dis¬ 
placements  w'  and  w  for  the  basic  state  that  are  non-zero. 

It  is  possible  to  obtain  a  complete  stiffness  matrix 
based  on  the  split  type  solution  in  the  following  manner.  Both 
the  displacements  and  stress  resultants  at  the  top  and  bottom 
edges  of  the  shell  element  can  be  computed  corresponding  to  the 
edge  effects  state  by  EQUATIONS  (4-18)  and  (4-19).  The  necessary 
coefficients  are  tabulated  in  TABLES  4.4  and  4.5.  In  a  similar 
manner  the  displacements  and  stress  resultants  at  the  edges 
corresponding  to  the  basic  state  can  be  obtained  by  EQUATIONS 
(4-24a)  and  ( 4- 24b )  and  the  required  coefficients  are  tabulated 
■fin  TABLES  4.7  and  4.8.  Using  the  values  of  the  stress  resultants 
thus  obtained  a  matrix  [SM]  can  be  written.  This  matrix  will  be 
8x8  with  the  first  four  columns  representing  the  stress  resultants 
for  the  displacements  applied  at  the  top  edge.  The  elements  are 
tabulated  in  TABLE  4.9.  The  elements  in  the  last  four  columns 
will  be  numerically  equal  to  the  corresponding  elements  in  the 
first  four  columns  with  appropriate  signs  in  accordance  with 
Flligge's  sign  convention.  It  may  be  noted  that  the  format  of 
of  [SM]  ,  that  is  the  labelling  of  the  rows  and  columns  resembles 
that  of  a  stiffness  matrix.  However,  although,  the  elements  of 
[  SM]  are  stress  resultants  at  the  edges,  [SM]  is  not  the  desired 
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stiffness  matrix  since  the  interaction  of  the  displacements 
caused  by  the  two  systems  has  not  been  considered.  To  trans¬ 
form  [sm]  to  the  desired  stiffness  matrix  it  is  necessary  to 
correct  the  elements  by  eliminating  the  non-zero  displacements  in 
each  system  caused  by  the  other  system.  This  is  done  as  follws. 

The  matrix  formed  from  the  elements  corresponding  to 
the  boundary  conditions  leading  to  the  stiffness  matrix  must  be  a 
unit  matrix.  In  a  manner  similar  to  that  adopted  for  writing 
a  matrix  [D]  corresponding  to  the  displacements  can  be  written 
using  EQUATIONS  (4-18) , (4-19) , (4-24a)  and  (4-24b)  and  these 
elements  are  tabulated  in  TABLE  4.10.  Any  column  in  this  matrix 
[D]  represents  the  applied  edge  displacements  and  the  resulting 
edge  displacements  considered  in  an  order  which  is  consistent  with 
that  of  [sm]  .  To  form  a  stiffness  matrix,  it  is  necessary  that  the 
matrix  [d]  be  converted  into  a  unit  matrix.  It  is  noted  that  all 
the  terms  on  the  major  diagonal  already  have  a  value  of  unity. 

The  displacements  in  the  first  column  of  the  matrix  [d! 
are  obtained  from  EQUATION  (4-18)  corresponding  to  the  edge  effects 
solution.  Therefore  the  only  non-zero  elements  other  than  on  the 
major  diagonal  will  be  displacements  corresponding  to  the  basic 
state.  These  elements  can  be  made  zero  by  applying  equal  and  opposite 
corrective  displacements.  Stress  resultants  corresponding  to  these 
displacements  can  be  computed  using  the  expressions  for  the  basic 
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state  given  in  ABLES  4.7  and  4.8.  The  values  are  added  to  the 
corresponding  elements  in  the  first  column  of  [ SM ]  .  This  proce¬ 
dure  is  repeated  for  the  second  column.  Correction  for  the  third 
and  fourth  columns  are  due  to  displacements  corresponding  to  the 
edge  effects  state.  The  corrections  are  obtained  using  TABLES  1 
4.4  and  4.5  in  a  similar  manner  as  for  the  first  column.  The 
corrections  applied  to  [sm]  due  to  the  membrane  displacements  in 
the  first  two  columns  of  [ D]  introduce  non-zero  bending  displace¬ 
ments  which  will  be  smaller  than  the  applied  displacements. 

Similarly  non-zero  membrane  displacements  will  be  introduced 
when  applying  the  corrections  due  to  the  bending  displacements  . 
in  the  third  and  fourth  columns  of  [d]  .  New  corrective  displace¬ 
ments  can  be  applied  to  balance  the  non  zero  displacements 
resulting  from  the  first  cycle  of  corrections.  This  procedure 
can  be  repeated  until  the  displacements  introduced  by  applying  the 
correction  are  to  any  desired  degree  of  smallness.  After  a  suffi¬ 
cient  number  of  iterative  cycles  the  elements  of  [sm]  will  approach 
those  of  the  desired  stiffness  matrix  [sm]  . 

Fortunately  it  is  not  necessary  to  use  this  iterative  proce¬ 
dure  to  obtain  a  stiffness  matrix  corresponding  to  the  split  type 
solutions,  since  explicit  solutions  can  be  written  for  all  elements 
of  the  stiffness  matrix  [SM]  in  terms  of  the  uncorrected  elements 
of  [sm]  and  [d]  .  These  expressions  are  given  in  TABLE  4.11  and 
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correspond  to  the  sign  convention  given  by  Flligge. 

The  expressions  given  in  TABLE  4.11  were  obtained  by 
applying  the  first  cycle  of  corrections  algebraical ly .  The  result¬ 
ing  elements  for  an  extreme  range  of  shell  parameters  were 
compared  with  the  corresponding  elements  of  the  stiffness  matrix 
obtained  by  Flugges's  exact  theory,  described  in  SECTION  4.2.  It 
was  noted  that  in  the  two  sub  matrices  corresponding  to  any  edge 
the  diagonal  terms,  terms  correspond! ng  to  w'  and  u  and  the  term  T 

A 

corresponding  to  w  agree  sufficiently  close  with  the  exact  values 
to  be  used  without  any  further • correction .  Since  the  final  matrix 
is  symmetrical  all  other  elements  can  be  written  from  symmetry.  A 
numerical  example  for  a  specific  shell  geometry  and  loading  harmonic 
is  given  as  Example  5.1  in  CHAPTER  V. 
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4 . 5  BASED  ON  THE  EIGHTH  ORDER  CHARACTERISTIC  EQUATION  FOR 

HIGH  VALUES  OF  m 

It  is  clear  from  the  discussion  in  CHAPTER  III  that  for 
high  values  of  m5the  general  solution  based  on  Donnell's  theory 
will  give  sufficiently  accurate  results.  Donnell's  eighth 
order  characteri sti c  equation  is  given  by  EQUATION (4-27) 

(A2  -  m2)4  +  1  A4  =  0  (4.27) 

This  equation  has  been  already  presented  in  CHAPTER  I  as 
EQUATION  (1-15).  For  high  values  of  m  the  solution  may  be 
simplified  by  ignoring  the  effect  of  the  far  edge.  The  solution 
for  w  can  therefore  be  written  in  the  form  given  by  EQUATION C3-5) . 
The  constants  C’-jto  C^are  evaluated  from  the  boundary  conditions 
at  one  edge. 

Because  the  number  of  constants  is  reduced  to  four  the 
general  procedure  explained  in  SECTION  4.2  for  the  exact  theory 
of  Flugge  can  be  used  to  establish  the  influence  coefficients. 

For  the  evaluation  of  the  constants  the  four  boundary  conditions 
correspondi ng  to  w' ,  w,  u  and  v  at  the  edge  y  =  o  are  considered. 

A  unit  value  of  each  of  these  displacements  is  applied  keeping 
the  remainder  of  the  displacements  as  zero.  Using  EQUATION  (3-1) 
any  quantity  7  either  a  displacement  or  a  stress  resultant  can 
be  written  as 

« 

f  =  c [a i C i  +  ^2^2!  ^  a3^3  ^ 


(4-28) 
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[C^  = 


'X1 

1 


a 


1 


yl 

0 

a  o 


-x2 

1 


a. 


u2 

0 

a  „ 


(4-30) 
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Using  this  equation,  the  matrix  equation  for  the  boundary 
conditions  can  be  written  in  the  following  form  which  has 
been  given  earlier  as  EQUATION  (4-1)  i 


[C,][C]  =  [I]  • 


(4-29) 


The  coefficient  matrix  [C-j]  is  given  by  EQUATION  (4-30).  Each 
column  vector  of  [C]  consists  of  the  constants  to  and 
each  column  vector  of  [I]  represents  the  displacement  boundary 
condition  for  each  case  of  the  applied  unit  edge  displacement. 
The  matrix  [C],  obtained  from  the  solution  of  EQUATIONS  (4-29) 
is  given  by  the  EQUATIONS  (4-31).  The  stress  displacement 
relations  on  which  Donnell's  theory  is  based  are  identical  to 
the  corresponding  simplified  expressions  of  Flligge.  Therefore 
TABLE  4.2  giving  the  values  of  the  characteristic  coefficients 
can  be  used  for  Donnell's  theory  also.  Using  this  table  and 
EQUATION  (4-28)  the  coefficient  matrix  [C^]  as  given  in 
EQUATION  (4-32)  for  the  stress  resultants  M  ,  S  ,  N  and  T 

XXX  X 

correspondi ng  to  the  displacements  w' ,  w,  u  and  v  is  obtained. 

The  stiffness  influence  coefficients  are  obtained  by  the 
following  equation. 


[K]  =  [C2][C] 


(4-33) 
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The  elements  in  the  stiffness  matrix  [K]  and  the  in¬ 
ternal  stress  resultants  in  the  shell  can  be  found  in  the 
following  manner.  The  elements  in  [Cl  and  [C^]  depend  only 
on  the  quantities  x  j>  y  •  =»  a.,  3*.  Explicit  expressions,  for 
these  derived  by  Hoff  based  on  Donnell's  theory,  have  been 
given  in  SECTIONS  A. 2  and  A. 6  of  APPENDIX  A.  Substituting 
these  values  into  EQUATIONS  (4-31)  and  (4-32),  [C]  and  [C2] 
can  be  obtained.  The  elements  in  [K]  are  found  by  premul tiplyi ng 
[C]  by  [C2].  The  internal  stress  resultants  for  any  case  of 
unit  applied  edge  displacement  can  be  found  using  EQUATION 
(4-28).  The  constants  C^  to  C^  in  this  equation  are  given  by 
the  particular  column  vector  of  [C]  correspondi ng  to  the  unit 
displacement  considered.  The  quantities  a-j  and  a2  have  been 
given  in  TABLE  4.2.  Values  of  a^  and  a^  can  be  obtained  from 
the  same  table  by  changing  the  subscripts  of  x  and  y  from  1 
to  2  and  those  of  a  and  3  from  1  and  2  to  3  and  4  respecti vely . 

It  may  be  noted  that  in  the  above  method  Donnell's  theory 
has  been  used  without  making  any  further  approximations. 

However  since  the  sign  convention  for  the  general  shell  theory 
has  been  used  in  evaluating  the  elements  of  the  stiffness  matrix, 
it  will  be  symmetric  except  for  the  signs.  A  truly  symmetric 
matrix  can  be  obtained  by  adopting  the  stiffness  matrix  sign 
convention.  (See  FIGURE  4.2). 
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4.6  BASED  ON  SPLIT  TYPE  SOLUTION  FOR  VERY  HIGH  VALUES  OF  m 

For  very  high  values  of  m  the  shell  behaves  as  if  it 
were  a  thin  plate.  As  a  result  the  membrane  stress  resultants 
N  and  T  and  the  correspond! ng  displacements  u  and  v  may  be 

X  X 

determined  by  the  equations  of  generalised  plane  stress  and 
the  bending  stress  resultants  M  and  S  and  their  corresponding 

X  X 

displacements  w'  and  w  by  the  equations  of  plate  bending.  These 
equations  can  be  obtained  from  the  three  simpler  differential 
equations  expressed  in  terms  of  displacements  corresponding 
to  the  equations  of  equilibrium  in  the  three  co-ordinate 
directions  as  presented  by  Flugge  (  p.  219,  ref.  14.). 

This  is  accomplished  for  plane  stress  by  omitting  terms  involv¬ 
ing  w  and  the  loading  in  the  first  two  equilibrium  equations 
and  for  plate  bending  by  retaining  only  the  terms  involving 
the  derivative  of  w  in  the  third  equilibrium  equation. 

The  shell  problem  is  thus  split  into  two  independent 
states  and  the  analysis  leading  to  the  establishment  of  influence 
coefficients  is  done  by  considering  these  two  states  separately. 

GENERALISED  PLANE  STRESS 

The  generalised  plane  stress  equations  are  given  as 


+  U")  +  C-j  ( U "  +  V")  =  0 

+  V")  +  c,  (u“  +  V'1  )  =  0 

(4-34) 

where  Ct  =  (1+u) 

1  TT^JT 
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Taking  the  solutions  for  u  and  v  as  u  cos  m<|>  and  v  sin  m<j> 
respectively ,  the  above  equations  can  be  reduced  to  the  following 
form. 


u"  ( 1  +  c.j )  -  ni  u  +  mc-|  v 1  =  o 


v' 


_  _ 

m  (1  +  c-j )  v  -  m  c i  u '  =  o 


(4-35a ,b) 


_  Ay  Ay 

The  solution  for  u  and  v  can  be  written  as  Ae  and  Be  respectively 
Substituting  these  in  the  above  equation  we  get  the  following 
characteristi c  equation. 


/,  2  2n 2 

(x  -  m  )  =  o 


(4-36) 


This  equation  has  four  roots  which  are  all  numerically  equal 
to  m,  two  of  these  being  repeated.  Therefore  the  solutions 
for  u  and  v  can  be  written  as, 


u  =  e"my[A1  +  A2y]  +  e,",J'[A3  +  A4y] 

7  =  e-my[B1  +  B2y]  +  e+my[B3  +  B4y] 

Since  m  is  very  large  the  second  term  in  these  solutions  can 
be  ignored  as  being  unbounded  as  y  increases.  Therefore, 

7  =  e~my[A,  +  A,y] 


trny. 


(4-37) 


1  2J 


=  e_my[B1  +  B2y] 


V 


(4-38) 
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Substituting  these  into  EQUATIONS  (4-35a)  , 

B-j  =  A]  -  (2  +  c-j  )  and 

mc-j 

B2  =  A2  •  (4-39 ) 

There  are  thus  two  arbitrary  constants  of  integration  A-j  and  A2 
which  are  to  be  evaluated  from  the  boundary  conditions.  The 
stiffness  influence  coefficients  for  the  generalised  plane  stress 
state  are  established  by  considering  the  boundary  conditions 
corresponding  to  u  and  v  as  follows, 


1 


1 


2+c 


1 


1 


(4-40) 


The  column  vector  £a.j  A2(j |  represents  the  constants  A-j  and  Ap 
for  the  boundary  conditions  u  =  1  and  v  =  o.  The  vector 
iAlv  A2v^  rePresents  the  same  constants  for  the  boundary  conditions 
u"  =  o  and  7=1.  EQUATIONS  (4-40)  are  solved  for  the  constants 
A-j  and  A2>  Using  EQUATIONS  (4-38)  and  (4-39)  the  membrane  dis¬ 
placements  corresponding  to  the  prescribed  sets  of  boundary 
conditions  can  be  evaluated. 
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The  membrane  stress  resultants  can  now  be  evaluated 
from  the  corresponding  displacements  u  and  v  obtained  from 
EQUATIONS  (4-38)  for  the  boundary  conditions  specified  as 
above  and  the  stress  displacement  relations  based  on  general¬ 
ised  plane  stress  theory.  These  relations  can  be  obtained  by 
omitting  the  terms  corresponding  to  the  bending  displacement 
in  the  stress  displacement  equations  as  given  by  FlUgge 
[EQUATIONS  12a-c,  p.  217,  ref.  14].  The  expressions  thus 
obtained  are  presented  in  TABLE  4-12. 

From  TABLE  4-12  the  values  of  the  stress  resultants 
at  y  =  o  can  be  found.  The  stiffness  matrix  ( S M 3 ]  correspondi ng 
to  the  membrane  displacements  u  and  v  can  finally  be  written  as, 


ms 

[SM3]= 


2 


L(l-u) 


(l-u). 


2 


(4-41) 


where  s  =  D(1 -u) 
r(3-u) 


PLATE  BENDING 

The  differential  equation  for  plate  bending  is 
wlv  +  2w*  1  * '  +  w"  =  0  •  (4-42) 

Taking  w  in  the  form  w  cos  1114  this  equation  is  reduced  to  the  form.? 


9  2- 

-  2m  w 


4- 

+  m  w  = 


(4-43) 


0 


81 


Introducing  w  =  C  e^,  we  get  from  the  above  equation 
the  following  characteristic  equation, 

(A2  -  m2)2  =  o  •  (4-44) 

This  equation  is  identically  the  same  as  EQUATION  (4-36) 
and  the  solution  for  w  is  of  the  form, 

w  =  e-^[Cl  +  C2y]  +  e+my[C3  +  C4y]  • 

Ignoring  the  unbounded  solution, 

w  =  e'my[C1  +  C2y]  •  (4-45) 

There  are  thus  two  arbitrary  constants  of  integration 
C-j  and  Cg  which  are  to  be  evaluated  from  the  boundary  conditions. 

The  stiffness  influence  coefficients  for  the  plate  bending  state 
are  established  by  considering  the  boundary  conditions  correspond¬ 
ing  to  w'  and  w.  The  following  equations  are  obtained  for  the 
constants  , 


-m 

f 

Clw ' 

Clw 

"l 

0 

_  1 

0 

-C2W' 

C2w- 

0 

1 

(4-46) 


The  column  vector  {C-|  ,  C^.j  represents  the  constants  C-j  and 
for  the  boundary  conditions  w'  ~  1  and  w  =  o.  The  vector 
l^lw  ^2w  }  rePresents  the  same  constants  for  the  boundary 
conditions  w'  =  o  and  w  =  1.  EQUATIONS  (4-46)  are  solved  for 
the  constants  C-j  and  C Using  EQUATIONS  (4-45)  the  bending  dis- 


[SM]  = 


-2m  K 
(l+u)m^i< 
0 
0 


(l+u)m2K 

2m3K 

0 

0 


0 

0 

2§m 

(1 -u)sm 


0 

0 

(1  -u)sm 
2sm 


M 


x 


x 


X 


T 


x 


(4-48) 
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placements  corresponding  to  the  prescribed  set  of  boundary  conditions 
can  be  evaluated. 


The  bending  stress  resultants  can  now  be  evaluated  from 
the  corresponding  displacement  w  obtained  from  EQUATION  (4-45)  for 
the.  boundary  conditions  specified  as  above  and  the  stress  dis¬ 
placement  relations  based  on  plate  bending  theory.  These  relations 
can  be  obtained  from  the  stress  displacement  equations  as  given 
by  Flu’gge  [  EQUATIONS  12  d-f,  p.  217,  ref.  14].  The  expressions 
thus  obtained  are  presented  in  TABLE  4.13. 


From  TABLE  4.13  the  values  of  the  stress  resultants  at 
y  =  o  can  be  found.  The  stiffness  matrix  [SM4]  corresponding 
to  the  bending  displacements  w '  and  w  can  finally  be  written  as, 


[SM4> 


w'  =  1 

w  =  1 

-2 

(l+u)m 

(l+u)m 

- 1 

CM 

E 

CM 

mK 


(4.47) 


COMPLETE  STIFFNESS  MATRIX 

Using  EQUATIONS  4-41  and  4-27  the  complete  stiffness 
matrix  considering  all  the  four  displacements  w',  w,  u  and  v 
can  be  written  as  given  in  EQUATION  4.48.  It  is  noted  that  the 
matrix  thus  obtained  is  truly  symmetric. 


TABLE  4.1 
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VALUES  OF  a-j  ,  «£>  3-j  >  3 ^ 


'/h 

k 

in 

al 

6i 

32 

0 

0.8333x1 0~3 

1 

Fliigge 

0.01240 

0.02765 

0.00110 

0.06334 

Donnel 1 

0.01163 

0.02797 

0.00538 

0.06335 

Simp! i fi ed 

0.01620 

0.02445 

0.00085 

0.06339 

3 

Fliigge 

0.01226 

0.06355 

0.01337 

0.18318 

Donnel 1 

0.00580 

0.06130 

0.01791 

0.18182 

Simp! i fi ed 

0.03164 

0.04297 

0.07820 

0.16035 

5 

Fliigge 

0.03095 

0.12528 

0.03875 

0.27316 

Donnel 1 

0.02162 

0.12285 

0.04307 

0.27316 

Simp! ified 

0.05766 

0.03435 

0.12199 

0.15090 

10 

Fliigge 

0.04831 

0.29466 

0.07583 

0.45888 

Donnel 1 

0.03582 

0.29160 

0.07626 

0.45762 

Simpl ified 

0.04855 

0.01888 

0.08822 

0.07451 

100 

0.8333xl0"5 

1 

Fliigge 

0.00615 

0.00663 

0.00001 

0.00634 

Donnel 1 

0.00626 

0.00652 

0.00000 

0.00635 

Simpl if i ed 

0.00613 

0.00665 

0.00005 

0.00634 

3 

Fliigge 

0.00512 

0.00769 

0.00017 

0.01902 

Donnel 1 

0.00524 

0.00759 

0.00023 

0.01899 

Simpl  i  fi  ed 

0.00487 

0.00777 

0.00104 

0.01899 

5 

Fliigge 

0.00325 

0.00994 

0.00095 

0.03158 

Donnel 1 

0.00338 

0.00985 

0.00104 

0.03157 

Simpl ified 

0.00216 

0.00969 

0.00458 

0.03123 

10 

Fliigge 

0.00259 

0.02090 

0.00620 

0.06021 

Donnel 1 

0.00237 

0.02082 

0.00634 

0.06015 

Simpl ified 

0.00931 

0.01288 

0.02478 

0.05157 

1000 

0.8333xl0"7 

1 

Fliigge 

0.00201 

0.00203 

0.00000 

0.00063 

Donnel 1 

0.00202 

0.00202 

0.00000 

0.00063 

Simpl ified 

0.00201 

0.00203 

0.00000 

0.00063 

3 

Fliigge 

0.00198 

0.00206 

0.00000 

0.00190 

Donnel 1 

0.00198 

0.00206 

0.00000 

0.00190 

Simpl ified 

0.00197 

0.00207 

0.00001 

0.00190 

5 

Fliigge 

0.00194 

0.00213 

0.00000 

0.00317 

Donnel 1 

0.00192 

0.00212 

0.00001 

0.00317 

Simpl ified 

0.00190 

0.00214 

0.00004 

0.00317 

10 

FI  iigge 

0.00161 

0.00245 

0.00008 

0.00634 

Donnel 1 

0.00162 

0.00245 

0.00008 

0.00634 

Simpl ified 

0.00152 

0.00247 

0.00037 

0.00632 

/  / 

Note  :  Simplified 
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APPROXIMATE  CHARACTERISTIC  COEFFICIENTS 
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TABLE  4.3 

VALUES  OF  ax  AND  a2  WITH  xi  =  HI 


f 

c 

ai 

a2 

T  factor 

w 

1 

1 

0 

cos 

w' 

1 

"X1 

X1 

cos 

u 

1 

al 

a2 

cos 

V 

1 

^1 

?2 

si  n 

T 

D 

1+U13-J  -ux-|  (a-|+a2) 

m32“uX-]  (a2"ai ) 

cos 

“x 

D 

u-x-j  (a-|  +«2  )+om3^ 

-x-j  (“2“ai  )+um^2 

cos 

V 

D(l-o)/2 

— mot  i  _x  i  ( 3 1 +  3  ^ ) 

-ma^-x-j  ( 32-3i  ) 

sin 

Tx 

D(l-o)/2 

-ma-j-x-]  (£-|+32) 

-ma^-x-]  ( 3 2 - 3 1  ) 

si  n 

T 

Mx 

K 

2 

-m 

»2oXl2 

cos 

K 

2 

-um 

-2o2 

cos 

‘Tx 

K(l-o)/2 

2mx-j 

-2mxi 

si  n 

T 

K 

m(m2-l ) 

2mx-|2 

cos 

T 

K/2 

2Xl(m2+2Xl2) 

- 2x i  (m2-2X-j 2 ) 

cos 

Sx 

K/2 

2Xl[(2-u)ni2+2Xl2] 

- 2x i  C  (2-u)m2-2x-| 2] 

cos 

Mote 

•  N 

•  XT 

Tx  ’  Hx*  = 

M 

Tx 
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TABLE  4.4 


CHARACTERISTIC  COEFFICIENTS 
FOR  w'  =  1  AND  w  =  o  AT  TOP 


f 

c 

si 

S2 

w 

1 

0 

1 

w1 

1 

1 

-1 

u" 

1 

a2/x-| 

T/xi 

V 

1 

32/X] 

Vxi 

\ 

D 

m32/X] _u^a2_al ^ 

( 1  +rn 3 -]  ) / X i  -o  (a-|  +a2 ) 

\ 

D 

-  (a^-a-j  )+um32/ X-j 

o(l+m3-|  )/xi  ~(ai+a2) 

V 

D(l-o)/2 

-ma2/xi -($2-3i  ) 

mot  -j  /  x  i  -  ( 3 1 +32 ) 

A. 

D(l-u)/2 

-ma2/x-]-(32-3-j  ) 

-ma i  / x -j  -  ( +32 ) 

2 , 

K 

-2ux-l 

-m  /x-| 

Mx 

K 

~2x1 

2, 

-um  /x-, 

M*. 

<bx 

K(l-o)/2 

-2m 

2m 

\ 

K 

2mx-, 

2 

m(m  -1 ) /x i 

K/2 

-2  (m2-2x-| 2) 

2  (m2+2x-j 2  )/x-j 

K/2 

-2[(2-u)m2-2Xl2] 

2[(2-u)m2+2X]2]/x1 

N 


Xcj) 


M 


<J>X 


Note: 
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TABLE  4.5 


CHARACTERISTIC  COEFFICIENTS 


FOR 

w1  =0  AMD  w  =  1 

AT  TOP 

f 

c 

S3 

S4 

¥ 

1 

1 

1 

w' 

1 

0 

-2X-i 

IT 

1 

a-j+o^ 

a-j  -ap 

v~ 

1 

V®2 

N, 

<t> 

D 

l+m($i  +  32)“2ox-]a2 

1  +m(  3-^  -32)_2vx-ja-] 

N 

X 

D 

u-2x-]a2+um(3-|  +  32) 

u-2x-ja-|+um(3i-32) 

D(l-u)/2 

-m(a-j+a2)"2x-| 

-ni(a-| -a^ )  “ 2 Xi  3-j 

T 

X 

D(l-u)/2 

-m(a-|+a2)-2x-|  ^2 

-m(ai -a2 ) ~ 2x i  3i 

. _ _ 

— 

2 

2 

M. 

4 

K 

-m  -2uXl 

-ill  +2oXl 

V 

K(l-u)/2 

0 

4mx-| 

% 

K 

m(m2-l  )+2mx-|  2 

m(m2-l  )-2inx-| 2 

2 

<4 

K/2 

8x,3 

4x-|m 

S 

X 

K/2 

V 

4X-,  (2-o)m2 

Mote: 

W  =  N*. 

X({)  (J)X 

5  M  ,  =  FT 
x<j>  Tx 
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TABLE  4.6 

VALUES  OF  THE  CONSTANTS  IN  THE 
BASIC  SOLUTION 


constants  Case  1  u  =1  ,v=o  at  top  Case  2  um  =0 ,  vm^ "1  at  top 

ITm=o ,7=o  at  bottom  um=^,vm=0  bottom 


A 

m 

1 

2G 

QG 

B 

-[m2Q2  +  2(1 +o)] 

-m 

3 

2G 

QG 

C 

1 

0 

D 

0 

1 

G  =  [m2Q2  +  2 ( 1 +u ) ] 

12 
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TABLE  4.7 

CHARACTERISTIC  COEFFICIENTS  FOR  7=1  AND  7=0, 
AT  TOP  AND  7=7=0  AT  BOTTOM 


f 

c 

F1 

F2 

F3 

F4 

w 

1 

-m4/12 

m2b/2Q 

m2[l+u/2-a] 

ub/Q 

w 1 

1 

0 

-m4/4 

m2b/Q 

m2[l+u/2-a] 

7 

1 

0 

m2/4 

-b/Q 

a 

7 

1 

m3/l  2 

-bm/2Q 

am-m(l+o) 

0 

\ 

D 

0 

0 

0 

2 

0 

Nx 

D 

0 

0 

m  /2 

-b/Q 

=  'VJ 

|  D(l-o)/2 

0 

0 

0 

-m/2 

,  b  =  m2Q2/3  +  2  (1+u) 


Note:  -  a  =  m2Q2/12  +  2  (1+u) 
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TABLE  4.8 

CHARACTERISTIC  COEFFICIENTS  FOR  u=o  AND  7=1 


AT  TOP  AND 

u=v=0 

AT  BOTTOM 

f 

c 

Fr 

5 

F6 

F7 

F8 

w 

1 

-m3/6Q 

m3/4 

(2+u)m/Q 

mu/ 2 -am 

w' 

1 

0 

-m3/2Q 

m3/2 

(2+U))m/Q 

u 

1 

0 

m/2Q 

-m/2 

0 

7 

1 

m3/6Q 

-in2/ 4 

-2(l+«)/Q 

a 

N. 

<P 

D 

0 

0 

0 

Q 

7 

X 

D 

0 

0 

m/Q 

-m/2 

V 
=  V 

D(l-. 

0/2  0 

0 

0 

-2/Q 

Note  :  a  =  m2Q2/12  +  2  (l+o)  ,  b  =  m2Q2/3  +  2  (1+v) 
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Note:  a  =  m  Q  /1 2  +  2(l+u)  ?  b  =  m  Q  /3  +  2  (1+u) 
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TABLE  4-10 

ELEMENTS  OF  THE  MATRIX  [D] 


D  ( i ,  j  ) 

Values  D(i , j ) 

D(1 ,1 ) 

1 

D(3,l) 

a2^xl 

0(4,1) 

VX1 

D (2,2) 

1 

0(3,2) 

CL  -j  “  CL 

0(4,2) 

V?2 

D(1 ,3) 

-[m2(l+u)/2-a]/a 

0(2,3) 

o  b/Q  a 

0(3,3) 

1 

D(1 ,4) 

-(2+\j)m/Q  a 

0(2,4) 

[om/2  -  a  m ]/  a 

0(4,4) 

1 

0(5,3) 

[m4Q2/4  m2  (b+l+o/2-b)]/b 

0(6,3) 

-m4Q3/l 2  +  [b/2+( 1 +u/2-a) ]  m2Q+u  b/Q 

0(5,4) 

-(2+u)m/Q  a 

D  ( 6 ,4 ) 

[m3Q2/l 2+(2+u)m  +  m/2  -  a  m]/  a 

Note:  Only  the  left  half  of  the  matrix  is  given. 
Values  not  shown  in  this  half  are  zero. 
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TABLE  4.11 

ELEMENTS  IN  THE  FINAL  MATRIX  [SM] 


SM(i,j)  Values  of  SM(i,j) 


SM ( 1 , 1  ) 

SM ( 1 ,2)=SM(2 ,1 ) 
SM ( 1  ,3)=SM(3,1 ) 
SM ( 1 ,4)=SM(4,1 ) 
SM(2,2) 

SM(2,3)=SM(3,2) 
SM(2 ,4)=SM(4,2) 
SM (3,3) 

SM(3,4)=SM(4>3) 
SM  ( 4 , 4 ) 

SM (5,1  ) 

SM (5,2)=  SM (6,1 ) 
SM(5,3)=SM(3,5) 
SM(5,4)=SM(8,1) 
SM(6,2) 

SM(6,3)=SM(7,1 ) 
SM(6,4)=SM(8,2) 
SM ( 7 , 3 ) 

SM(7,4)=SM(8,3) 
SM  ( 8 , 4 ) 


SM  (1,1) 
SM ( 2 , 1 ) 
SM ( 1 ,3) 
SM (4,1) 
SM (2,2) 
SM  ( 2 , 3 ) 
SM  ( 4 , 2 ) 
SM(3,3) 
SM  ( 3 , 4 ) 
SM(4,4) 
0 
0 

SM  ( 5 , 3 ) 
SM(8, 1 ) 
0 

SM ( 6 , 3 ) 
SM(8,2) 
SM ( 7 , 3 ) 
SM(8,3) 
SM(8,4) 


D (3 ,1 ).  SM ( 1  ,3) 
D ( 3 , 1 ) .  SM  ( 2 , 3 ) 
D  ( 1 » 3 ) .  SM  ( 1  , 1 ) 
D  ( 3 , 1 ) .  SM  ( 4 , 3 ) 
D ( 3 , 2 ) - SM (2,3) 
D(1 ,3)  -  SM  (2,1 ) 
D(3,2).  SM(4,3) 
D  ( 1  , 3 )  •  SM  ( 3 , 1 ) 
D  ( 1  ,3)  *  SM  (4,1 ) 
D(1,4).SM(4J) 


D(5,3).  SM  ( 1 ,1 ) 
D  ( 3 , 1 )  •  SM  ( 8 , 3 ) 


D(4,l) 

SM ( 1 ,4) 

D(4,l) 

SM(2,4) 

D (2,3) 

SM ( 1 ,2) 

D(4,l) 

SM (4,4) 

D(4,2) 

SM  ( 2 , 4 ) 

D(2,3) 

SM(2 ,2 ) 

D(4,2) 

SM (4 ,4) 

□(2,3) 

SM (3,2) 

D(2,3) 

SM (4,2) 

D(2,4) 

SM  ( 4 , 2 ) 

D(6,3). 

SM ( 1 ,2) 

D(4,l) 

SM  ( 8 , 4 ) 

D (5,3).  SM(2,1 )  +  D ( 6 , 3 )  •  SM  ( 2 , 2 ) 
D(3,2).SM(8>3)  -  D  ( 4 , 2 )  •  SM  ( 8 , 4 ) 
D  ( 5 , 3 ) .  SM  ( 3 , 1 )  -  D  ( 6 , 3 ) .  SM  ( 3 , 2 ) 
D (5 ,3) • SM (4 , 1 )  +  D(6,3)*SM(4,2) 
D (5,4) • SM (4,1 )  +  D(6,4).SM(4,2) 
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TABLE  4.12 

MEMBRANE  STRESS  RESULTANTS 
FOR  VERY  HIGH  VALUES  OF  m 


f 

c 

Case  1 :  -  u=l ,  v=0 

Case  2:  u=0,  v=l 

"x 

D 

[(c-m)-mc(l-u)y]g 

[-c+um+mc(l-u)y]g 

% 

D 

[u(c*-m)+mc"(l  -u)y]g 

[  (m-uc’)-nic’(l  -u)y]g 

NxM 

D(l- 

-u)/2  [(c-m)-2mc  y]g 

[  -  (c+rn)+2mc  y]g 

V! 

r 

Note 

:  c 

=  m(l+u)/ (3-u)  ,  g  =  e“my 
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TABLE  4.13 

BENDING  STRESS  RESULTANTS 
FOR  VERY  HIGH  VALUES  OF  m 


f 

c 

Case  1:  -  w 1 =1  w=0 

Case  2:  -  w'=0,  w=l 

Mx 

K 

[-2m+m2(l-u)y]  g 

[- (l+u)+m(l -u)y]  m2g 

\ 

K 

[l-m(l-u)y]m  g 

2m3g 

K 

[-2\jm-m  (l-u)y]g 

[-  (l+u)-m(l  -u)y]m2g 

% 

K 

2m  2  g 

[(l+u)+m(l-u)+m(l-u)  (l-m)y]f 

M , 

Tx 

K 

[- (1 -u)+(l -u)my]mg 

(l-u)m  yg 

Note 

-my 
g  =  c  J 

TOP 
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CHAPTER  V 


CRITERIA  FOR  THE  PROPOSED  APPROXIMATE  METHODS  AND 

NUMERICAL  EXAMPLES. 

5.1  INTRODUCTION 

In  this  chapter  criteria  for  the  range  of  applicability 
of  the  approximate  methods  for  establishing  the  stiffness  in¬ 
fluence  coefficients  described  in  CHAPTER  IV  are  presented.  The 
criteria  are  based  not  only  on  the  values  of  the  loading  harmonic 
m,  but  also  on  the  shell  geometric  parameters  k  and  Q.  For  ease 
of  reference  the  different  approximate  procedures  are  classified 
as  follows, 

f/ethod  1  :  solution  based  on  splitting  of  the  differential 
equation  for  low  values  of  m  (SECTION  4.4). 

Method  2  :  solution  based  on  simplified  Donnell's  theory 
for  high  values  of  m  (SECTION  4.5). 

Method  3  :  solution  based  on  equations  of  generalised 
plane  stress  and  plate  bending  for  very  high  values  of  m.  (SECTION  4.6). 

Examples  are  included  to  illustrate  the  regions  of 
applicability  of  the  above  methods.  Example  1  considers  the  analysis 
of  a  steel  bin  for  wind  effects  which  involves  low  values  of  m 
and  Example  2  considers  the  analysis  of  a  water  tank  on  isolated 
supports  which  involves  high  and  very  high  values  of  m.  The  range 
of  applicability  of  Method  3  is  discussed  in  this  latter  example. 
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5.2 


RANGE  OF  APPLICABILITY  OF  METHODS  1  AND  2 
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The  criteria  defining  the  range  of  applicability  of 
Methods  1  and  2  were  established  by  comparing  the  final  stress 
resultants  obtained  by  these  methods  with  those  obtained  using 
the  exact  theory  proposed  by  Flugge.  For  this  purpose  a  closed 
cylindrical  shell  free  at  the  top  and  fixed  at  the  bottom  was 
considered.  The  shell  was  analysed  for  wind  effects  in  which 
the  loading  was  taken  in  the  form  of  a  Fourier  series  as  given 
by  the  following  equation, 

q  =  qQ  +  q-jCOScj)  +  q2COs2<j>  +  q^cosSf  +  ... 

For  each  harmonic  beyond  the  first  term  of  this  series  the 
analysis  was  obtained  by  Flugge's  exact  theory,  Method  1  with 
one  cycle  correction  and  Method  2  neglecting  the  far  edge  effect. 

The  effects  of  variations  in  k  and  Q  on  the  final  stress 
resultants  M  ,  S  ,  N  and  T  at  the  support  are  presented  graph- 
ically.  The  comparison  of  the  results  obtained  by  the  exact 
theory  and  Method  1  are  given  in  FIGURES  5.1  to  5.3  and  values 
obtained  by  the  exact  theory  and  method  2  are  given  in  FIGURES  5.4 
to  5.6. 

The  figures  indicate  the  following.  For  a  given  shell, 
as  m  increases  Method  1  becomes  less  accurate  and  Method  2 
becomes  more  accurate.  However  when  either  the  thickness  or  the 
length  of  the  shell  is  reduced  (k  or  Q  is  reduced)  Methdd  2 
becomes  less  accurate.  The  reasons  are  as  follows.  From  SECTION  4.4 
it  is  clear  that  the  splitting  of  the  original  eighth  order 
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characteristi c  equation  for  low  values  of  m  is  more  justified 
when  k  is  very  small.  Also  in  Method  1  the  basic  solution  is 
replaced  by  the  membrane  solution.  It  can  be  seen  from  the 
inequality  given  by  EQUATION  (4-22)  that  when  values  of  either 
k  or  Q  become  smaller,  the  value  of  m  for  which  this  inequality 
is  satisfied  becomes  greater.  Therefore  for  low  values  of  k  or 
Q  the  assumptions  on  which  Method  1  is  based  become  more  just¬ 
ified  and  the  accuracy  increases.  However  Method  2  does  not 
give  good  results  for  these  conditions  as  the  assumptions  on 
which  this  method  is  based  do  not  permit  an  accurate  evaluation 
of  x£  and  y^.  The  value  of  x^  1S  very  sma  1 1  and  cause  slowly 
damped  solutions.  For  any  given  value  of  m,  when  k  is  large, 
also  becomes  large  and  is  evaluated  more  accurately.  In  addition1, 
when  Q  is  large,  the  damping  of  the  x^"  solution  is  faster  for 
any  given  value  of  m  and  k  and  the  assumption  of  neglecting  the 
far  edge  becomes  more  justified.  Therefore  for  high  values  of  k 
and  Q  the  accuracy  of  Method  2  increases. 

Based  on  the  above  discussion  regarding  the  final  stress 
resultants  for  any  loading  harmonic,  it  is  seen  that  the  criteria 
in  terms  of  m  defining  the  ranges  of  validity  of  Methods  1  and  2 
are  functions  of  not  only  m  but  also  k  and  Q.  These  are  therefore 
established  on  the  basis  of  the  maximum  tolerance  that  can  be 
permitted  in  the  final  stress  resultants  for  any  values  of  k  and  Q. 

From  FIGURES  5.1  to  5.3,  the  values  given  by  Method  1 
become  much  greater  than  the  exact  values  for  higher  values  of  m 
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and  the  maximum  percentage  difference  is  observed  in  T  .  There- 

X 

fore  the  maximum  value  of  in  for  which  Method  1  is  valid  is  found 
by  specifying  a  maximum  tolerance  for  the  difference  in  T  . 

X 

The  tolerance  arbitrarily  specified  was  30  percent  for  the 
worst  term  of  T  .  It  should  be  noted  that  this  difference  of 

X 

30  percent  specified  for  T  corresponds  to  the  difference  for 

A 

only  a  particular  loading  harmonic  and  not  the  difference  for 
the  complete  analysis  for  wind  load  which  is  expressed  as  a 
trigonometric  polynomial.  The  following  equation  has  been  proposed 
by  Soare  (21)  for  wind  load  on  a  smooth  surface. 

w  =  g  [-0.804  +  0.140cos.4+  1.380cos2<f> 
w 

+  0.490cos3(j)  -  0.318cos4(f)]  (5-1) 

Where  g^  is  the  basic  dynamic  pressure.  In  this  equation  the 
coefficient  is  maximum  for  the  term  for  which  m  =  2  and  the 
coefficients  for  higher  loading  harmonics  become  increasingly 
smaller.  Because  of  this  in  a  complete  analysis  for  wind  load 
the  net  difference  is  expected  to  be  much  smaller  than  the 
maximum  error  in  any  given  term.  This  is  illustrated  later  in 
EXAMPLE  1,  SECTION  5.3  where  the  maximum  difference  of  30  percent 
occurs  for  the;,term  T  for  which  m  =  3,  but  the  difference 

X 

from  summing  the  terms  up  to  and  including  this  term  is  only 
18  percent.  The  difference  for  other  stress  resultants  will 
of  course  be  smaller. 
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Thus  the  maximum  tolerance  of  30  percent  in  a  given 
term  of  T  is  therefore  reasonable  and  can  be  accepted.  On 
this  basis,  considering  the  stress  resultant  T  the  maximum  value 

X 

of  m  for  which  Method  1  is  valid  is  found  from  FIGURES  5.1  to  5.3. 
With  these  maximum  values  FIGURE  5.7  is  obtained  which  shows  the 
relation  among  the  maximum  value  of  m,  k  and  Q, 

The  criterion  for  Method  2  is  established  as  follows. 
Referring  to  FIGURES  5.4  to  5.6  it  can  be  seen  that  the  stress 
resultants  given  by  this  method  are  always  less  than  the  exact 
values  for  low  values  of  m,  the  maximum  difference  being  in  N  . 

X 

As  the  values  were  not  conservative  the  minimum  value  of  m 
which  gave  conservative  values  of  the  worst  term  for  the  stress 
resultant  N  was  noted  from  FIGURES  5.4  to  5.6  for  each  value 

X 

of  k  and  Q.  FIGURE  5.8  thus  obtained  shows  the  relation  among 
the  minimum  value  of  m,  k  and  Q.  A  similar  figure  could  be 
obtained  by  permitting  non-conservative  value  of  N  ,  say, corres- 

A 

ponding  to  a  difference  of  10  percent  in  the  worst  term  N  . 

A 

This  would  have  the  effect  of  changing  the  criterion  increasing 
the  range  of  applicability  of  Method  2. 

The  criteria  for  the  ranges  of  validity  of  Methods  1  and 
2  are  thus  given  by  FIGURES  5.7  and  5.8.  For  any  given  values 
of  k  and  Q  the  maximum  and  minimum  values  of  m  for  which,  re se 
pectively,  Methods  1  and  2  are  valid  can  be  found  using  these 
figures.  The  use  of  these  figures  is  illustrated  in  EXAMPLE  1 
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given  in  SECTION  5.3. 

It  is  noted  from  FIGURES  5.4  to  5.6  that  Method  2  gives 
accurate  values  as  m  becomes  large.  Therefore  for  high 
values  of  m  Method  2  can  be  expected  to  give  results  which  are 
close  to  those  that  can  be  obtained  by  the  exact  theory. 

For  most  shells  of  dimensions  encountered  in  practice 
all  of  the  terms  in  the  series  for  wind  load  can  be  evaluated 
within  the  tolerance  established  by  the  above  criteria  by  either 
Method  1  or  by  Methods  1  and  2.  For  certain  cases  there  may 
be  a  particular  loading  harmonic  that  does  not  fall  within  the 
range  of  applicability  of  these  methods  as  established  by  the 
above  proposed  criteria.  In  these  cases  Method  1  should  be 
used.  This  means  that  the  maximum  difference  for  this  term  only 
will  exceed  the  maximum  difference  specified  for  a  term  on  which 
the  criterion  was  based,  but  the  result  will  always  be  on  the 
conservative  side.  The  final  values  of  the  stress  resultants 
will  be  sufficiently  accurate  for  design  purposes. 

5.3  NUMERICAL  EXAMPLES 

EXAMPLE  1. 

This  example  illustrates  the  use  of  Method  1  for  the 
analysis  of  a  vertical  cylindrical  steel  bin  for  wind  effects. 

The  structure  is  assumed  to  be  fully  fixed  at  the  bottom  and  free 
at  the  top.  Therefore  the  boundary  conditions  are  w'  =w=u=v=o 
at  the  bottom  and  Mv  =  SY  =  N  =  T  =  o  at  the  top.  The  design 

X  X  X  X 
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data  are 

Radius  r  =  25  ft 

Thickness  h  =  ^/Z  in 

Length  Z  =  125  ft 

Poisson's  ratio  u  =  0.3 

The  wind  load  is  taken  in  the  form  of  a  Fourier  series 
along  the  directrix  as  given  by  EQUATION  (5-1).  The  particular 
solution  for  the  surface  loading  due  to  wind  may  be  obtained 
by  considering  sufficient  terms  of  the  load  expressed  as  a 
Fourier  series  along  the  generator  of  the  shell.  However  in 
this  example  as  is  usual  in  the  analysis  of  shell  structures 
the  membrane  solution  has  been  taken  as  a  close  approximation 
for  the  particular  solution.  The  homogeneous  solution  which 
accounts  for  the  edge  effects  is  determined  from  the  stiffness 
matrix  established  by  Method  1. 

A  complete  analysis  for  the  stresses  in  the  structure 
requires  that  each  term  of  the  Fourier  series  for  the  loading 
along  the  directrix  be  evaluated  separately.  The  first  step 
is  to  determine  the  number  of  terms  in  this  series  for  which 
Method  1  is  applicable.  This  is  done  as  follows, 

P  =  -  =  0.0016667 
r 

Q  =  -  =  5 

x  r 

2 

k  =—  =  0.231 48x1 0-6 


107 


-0 . 2669E-07 

"  0 . 8411E-06 

0 . 3337E-06 

0 . 5297E-04 

0 . 1335E-06 

-0 . 4196E-05 

- 

0 . 1690E-05 

-0 . 1059E-03 

- 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

_ 

•  0  0.0 

•°  0.0 

. 1096E-02  -0 . 1524E-03 

.  1524E-03  -0 . 3049E-04 


.0  0.0  | 

.0  0.0  J 

. 4289E-03  0 . 1524E-03  | 

. 1524E-03  -0.3049E-04! 

(5-2) 


0 

0 

■0 

0 

0 

0 

0 

0 
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0. 100E01 

0.0 

0.3579E  01 

-0.8415E-0ll 

0.0 

0.1000E  01 

0.1972E  00 

-0.1973E  01 

-0 . 1532E-03 

0 . 9528E-02 

0.1000E  01 

0.0 

-0. 7368E-04 

0. 2309E-02 

0.0 

0.1000E  01 

0.0 

0.0 

-0.4207E  00 

-0 . 8415E-01 

0.0 

0.0 

-0 . 7721E-01 

-0 . 2744E-01 

0.0 

0.0 

O'.O 

0.0 

0.0 

0.0 

0.0 

0.0 

(5-3) 
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Referring  to  FIGURE  5.7,  for  these  values  of  k  and  Q, 

m  =3 
max 

This  means  that  Method  1  is  valid  for  evaluating  all  stress 
resultants  for  which  m  is  equal  to  3  or  less  within  the  criterion 
used  to  obtain  FIGURE  5.7.  This  can  also  be  demonstrated  by 
checking  the  value  of  p  given  by  EQUATION  (4-21)  and  reproduced 
here  for  reference, 

p  =  m^(m^-l )  P^Q^  <<  1 
12(47) 

Si  nee, with  the  values  of  P  and  Q  in  this  example,  p  is  equal  to 
0.0039  which  is  much  smaller  than  l,the  assumptions  on  which 
Method  1  is  based, are  valid  and  Method  1  can  be  used. 

To  illustrate  the  use  of  the  stiffness  matrix,  detailed 
calculations  are  shown  for  the  term  m  =  2  only.  The  terms  for 
m  =  1  and  3  are  evaluated  in  a  similar  manner. 

The  stiffness  matrix  for  m  =  2  is  established  using  the 
procedure  explained  in  SECTIONS  4.4.  From  EQUATIONS  (A-3)  and 
(A-5)  values  of  x-j  and  «q  *  c^,  3^  and  are  evaluated.  Sub¬ 
stituting  these  values  in  the  expressions  given  in  TABLES  4.9 
and  4.10  the  matrix  [SM]  and  the  matrix  of  unbalanced  displacements  ft)"] 


0 . 2670E-07 

0 . 8387E-06 

-0 . 7132E-07 

0 . 1665E-05 

0 . 8387E-06 

0 . 5297E-04 

-0 . 7433E-05 

0 . 1043E-03 

-0. 7132E-07 

-0 . 7433E-05 

0 . 1097E-02 

0 . 1378E-03 

0.1665E-05 

0 . 1043E-03 

0 . 1378E-03 

0 . 2362E-03 

0.0 

0.0 

-0. 7657E-07 

0 . 2559E-07 

0.0 

0.0 

-0 . 4443E-05 

0 . 1523E-05 

-0 . 7657E-07 

-0 . 4444E-05 

0 . 4294E-03 

0 . 1437E-03 

-  0 . 2559E-07 

-0 . 1523E-05 

0 . 1437E-03 

0 . 2749E-04 

0.2664E-07 
0.8378E-06 
-0.7125E-07 
0 . 1649E-05 

0 . 1348E-10 
0.7858E-09 
-0.7649E-07 
-0.2386E-07 


0 .8378E-06 
0.5287E-04 
-0 . 7424E-05 
0 . 1042E-03 

0 . 7858E-09 
0.4573E-07 
-0.4440E-05 
-0.1431E-05 


-0 .7124E-07 
-0.7424E-05 
0.1097E-02 
0 . 1380E-03 

-0.7649E-07 
-0 . 4440E-05 
0.4300E-03 
0.1439E-03 


0 . 1649E-05 
0 . 1042E-03 
0 . 1380E-03 
0 . 2360E-03 

-0 . 2386E-07 
-0.1431E-05 
0 . 1439E-03 
0.2466E-04 
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are  evaluated.  These  values  are  given  in  EQUATION  >  (5-2)  and 
(5-3).  The  elements  in  these  matrices  are  in  accordance  with 
the  sign  convention  used  in  the  general  shell  theory  and  as 
given  by  FIGURE  4.1.  Only  half  the  portion  of  the  matrices  [SM] 
and  [D]  are  given  as  the  correction  to  be  applied  to  [SM]  involve 
the  elements  in  these  halves  only.  The  correction  is  applied 
as  described  in  SECTION  (4.4.3).  The  corrected  matrix  is 
written  using  the  sign  convention  shown  in  EQUATION  (5-4).  Again 
only  half  of  the  final  stiffness  matrix  [SM]  is  shown  in 
EQUATION  (5-4)  as  the  other  half  can  be  written  using  EQUATION  (4-8). 
It  is  seen  that  [SM]  agrees  very  well  with  the  stiffness  matrix 
[Kg]  given  in  EQUATION  (5-5),  obtained  by  the  exact  theory. 

Once  the  stiffness  matrix  [SM]  is  established  for  a 
particular  harmonic  m  there  are  two  procedures  as  described  in 
CHAPTER  2  to  obtain  a  final  solution.  The  first  procedure  is  to 
use  the  classical  procedure  by  writing  the  homogeneous  solution 
and  the  particular  solution  and  superimposing  these.  The  unknown 
displacements  in  the  homogeneous  solution  are  evaluated  by  making 
the  general  solution  satisfy  the  specified  boundary  conditions. 

In  this  problem  there  are  eight  boundary  conditions,  four  expressed 
in  terms  of  force  at  the  top  edge  and  the  remaining  four  expressed 
in  terms  of  displacements  at  the  bottom  edge.  The  four  homo¬ 
geneous  displacements  at  the  bottom  edge  must  be  equal  and  opposite 
to  the  displacements  given  by  the  particular  solution  to  satisfy 
the  displacement  boundary  conditions  at  the  bottom  edge.  Therefore 


these  four  quantities  can  be  evaluated  directly.  There  are  now 
only  four  unknown  homogeneous  displacements  at  the  top  and  these 
are  evaluated  by  inverting  a  4x4  matrix  written  by  using  the 
upper  half  portion  of  the  stiffness  matrix  [SM]  and  the  four 
force  boundary  conditions  at  the  top. 

If  one  were  to  use  the  alternative  stiffness  matrix 
procedure  explained  in  CHAPTER  11  the  initial  forces  [x 1 J  as 
given  by  EQUATION  (2-12)  must  be  found.  In  this  equation  the 
stiffness  matrix  [SM]  must  be  substituted  for  [K].  The  quantities 
£d'p]  andpCp]  are  the  edge  displacements  and  forces  obtained 
from  the  particular  solution  or  the  membrane  solution.  Using 
EQUATION  (2-13)  the  unknown  displacements  given  by  the  vector 
[d1]  is  found.  The  elements  in  the  vector  [X'J  are  the  applied 
forces,  four  of  which  are  known  to  be  zero  at  the  top.  In  the 
vector  [d1]  the  four  elements  at  the  bottom  are  known  to  be 
zero  and  thus  there  are  only  four  unknowns  to  be  found  by  in¬ 
verting  a  4x4  matrix  as  before.  The  advantage  in  using  this 
method  is  that  it  permits  the  evaluation  of  the  net  displacements 
at  the  edges  directly. 

However  in  this  example  the  complete  analysis  was  done 
using  the  former  classical  approach.  TABLE  5.1  gives  the  values 
of  the  stress  resultants  at  the  bottom  of  the  shell  both  by 
Method  1  and  the  exact  theory.  The  analysis  for  the  fourth 
harmonic  was  also  done  by  Method  1  even  though  this  harmonic 
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was  outside  the  range  of  validity  established  for  Method  1. 

The  results  given  are  for  unit  value  of  the  basic  dynamic  pressure 
gw  (see  EQUATION  (5-1)).  It  is  seen  from  this  table  that 
even  though  the  error  in  T  is  30  percent  corresponding  to  the 

A 

value  of  m  =  3,  in  the  complete  analysis  the  error  in  T  is 
only  18  percent  and  further  this  is  on  the  conservative  side. 

The  errors  in  the  other  stress  resultants  are  less.  This  means 
that  Method  1  can  be  used  for  harmonics  beyond  the  maximum 
value  of  m  obtained  using  the  criterion  for  the  range  of  validity 
for  this  method.  However  the  magnitude  of  error  in  the  final 
stress  resultants  depends  on  the  relative  magnitudes  of  the  Fourier 
coefficients  in  the  series  expressing  the  load. 

EXAMPLE  2 


To  illustrate  the  application  of  Methods  2  and  3  a 
circular  tank  supported  on  6  columns  is  analysed.  It  is  assumed 
that  the  tank  has  been  first  analysed  by  the  ordinary  theory 
of  symmetrical ly  loaded  tanks  on  the  assumption  of  a  continuous 
support  along  the  circumference.  The  scope  of  the  calculation 
shown  below  is  to  find  the  additional  stresses  caused  by  the 
discontinuity  due  to  the  supports  being  isolated. 


The  design  data  are: 
Length  % 

Thickness  h 

Radius  r 

Column  width 


39  ft 

5.9  ins 

32.8  ft 

frr  =  4.3  ft 
24 


Column  distance 


Trr  =  34.3  ft 
3 

0.25 


115 


Poisson's  ratio  u 


Nx  from  the  symmetrical  load  analysis  is  distributed 
uniformly  along  the  bottom  edge.  The  total  load  is  33.6  nr  kips/ft. 
A  self  equilibrating  system  of  N  distributed  as  shown  in 

X 

FIGURE  5.9  is  applied  to  the  bottom  edge  to  eliminate  the  value 
of  N  between  columns.  The  load  on  the  column  is  assumed  to  be 
constant  over  the  column  width.  This  edge  load  N  is  assumed 
as  the  Fourier  series 


IT 

where  m  =  6n. 


sin  nrr_  cos  m<j> 
8 


The  other  boundary  conditions  are: 

v  =  w  =  w'  =  o. 

0  0  0 

The  influence  from  the  upper  edge  is  neglected. 

The  stiffness  matrix  can  be  evaluated  for  each  harmonic 
using  either  Method  2  or  3  using  the  procedure  explained  in 
SECTIONS  (4-5)  and  (4-6).  For  illustration,  EQUATION  (5-6a)  and 
( 5 -6b )  give  the  stiffness  matrices  obtained  by  Method  2  and 
Method  3  respectively  for  m  =  6. EQUATION  (5-7a)  and  (5-7b)  give 
the  corresponding  quantities  for  m  =  120.  It  is  seen  that  for 
the  very  high  value  of  m  (m  =  120)  the  bending  terms  correspond¬ 
ing  to  the  applied  bending  displacements  and  the  membrane  terms 
correspond!' ng  to  the  applied  membrane  displacements  obtained  by 
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[K]  = 


[K]  = 


0 . 571 3E-05 
0. 5743E-04 
0 . 7657E-04 
0 . 2423E-03 


0.5143E-04 
0.121 5E-02 
0. 7308E-03 
0 . 651 6E-02 


0 . 7657E-04 
0 . 7308E-01 
0 . 41 82E-01 
0.1315E-01 


0.2423E-01 
0 . 1 31 5E-01 
0.131 5E-01 
0.4185E-01 


0.8600E-05 

0 . 1 350E-04 

0.0 

0. 1350E-04 

0.1 296E-03 

0.0 

0.0 

0.0 

0 . 5236E-01 

0.0 

0.0 

0 . 1 963E-01 

(5-6a) 

0.0 

0.0 

0 . 1 964  E-01 
0.5236  E-01 


(5-6b) 
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[K]  = 


[K]  = 


"  0. 7200E-04 

0. 5400E-02 

0.6818E-05 

0 . 1 81 8E-04 

0 . 5400E-02 

0 . 1 037E-01 

0.1091 E-02 

0 . 7909E-02 

0 . 681 8E-05 

0.1091E-02 

0.1047E  01 

0 . 3927E-00 

-0.881 7E-05 

0. 7909E-02 

0.3927E  00 

0.1047E  01 

(5-7a) 


~ 0 . 7200E-04 

0 . 5400E-02 

0.0 

0.0 

0 . 5400E-02 

0 . 1 037E-01 

0.0 

0.0 

0.0 

0.0 

0.1047E 

01 

0.3927E-00 

0.0 

0.0 

0.3927E 

00 

0.1047E  01 

( 5  -  7b ) 


Method  3  agree  very  well  with  those  by  Method  2.  The  agreement 
is  not  good  for  the  lower  value  of  m  =  6. 


The  shell  is  fully  analysed  considering  the  first  20 

terms  to  study  the  nature  of  convergence  of  the  final  stress 

resultants.  Because  of  the  uncoupled  nature  of  the  stiffness 

matrix  in  Method  3  only  the  boundary  conditions  correspond!' ng  to 

Nxq  and  vq  are  to  be  considered  when  using  this  method  for  this 

problem  and  as  such  this  method  gives  values  of  only  Nx, 

and  Na. 

4> 

TABLE  5.2  gives  the  values  of  the  stress  resultants  at 
the  support  after  considering  20  terms.  It  is  seen  that  the 
values  of  the  membrane  stress  resultants  at  the  support  obtained 
by  Method  3  agree  exactly  with  the  results  of  Method  2. 


FIGURES  5.10  and  5.11  illustrate  the  convergence  of  the 


stress  resultants  M  and  N  at  the  support  and  at  a  distance  of 

X  X 

it/8  from  the  support.  The  values  were  obtained  by  Method  2. 

The  convergence  of  M  is  very  rapid  and  the  value  of  M  can  be 
obtained  by  considering  the  first  3  or  4  terms,  but  the  con¬ 
vergence  of  Nx  is  slow  and  that  of  very  slow  so  that  more 
than  20  terms  are  required.  At  &/8  from  the  support  however 
the  convergence  of  both  N  and  N  also  are  very  fast  and  only 

X  X(J) 

4  to  5  terms  are  required. 


FIGURE  5.12  shows  the  variation  of  the  stress  resultants 
along  the  length  of  the  shell,  these  values  being  obtained  by 
considering  20  terms  in  the  series.  It  is  seen  that  Mx  and  Nx 


- 
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become  damped  at  a  faster  rate  than  NX(j)  .  The  maximum  value  of 
NX(j)  occurs  at  a  distance  of  V6from  the  support  at  which  location 
fortunately  as  indicated  by  FIGURE  5.11  the  convergence  of 
is  also  rapid.  Therefore  this  maximum  value  can  be  obtained 
by  considering  only  a  few  terms  in  the  Fourier  series. 

Even  though  Donnell's  theory  does  not  g’ve  correct  values 
of  y-^and  y£  (  refer  to  page  14)  for  very  high  values  of  m,  Method  2 
which  is  based  on  Donnell's  theory  was  found  to  give  final  results 
within  an  accuracy  of  1  to  2  percent  when  compared  with  the  exact 
results  obtained  by  Flugge's  theory  .  This  indicates  that  Donnell's 
theory  can  be  used  in  place  of  the  exact  theory  for  problems  of 
this  nature.  From  Figure  5.12  it  is  seen  that  the  membrane  stress 
resultants  obtained  by  Method  3  agree  well  with  the  results  obtained 
by  Method  2.  Therefore  similar  problems  can  be  completely  analysed 
by  Method  3  for  the  membrane  stress  resultants.  The  advantage  of 
this  method  is  that  only  a  2x2  matrix  must  be  considered.  The 
bending  stress  resultants  are  obtained  by  Method  2  considering  the 
first  3  or  4  terms  in  the  series  depending  on  the  desired  convergence. 


120 


TABLE  5.1 

COMPARISON  OF  STRESS  RESULTANTS 
BY  EXACT  THEORY  AND  METHOD  1 . 


Exact  Theory 

Method 

1 

m 

Sx 

Nx 

Tx 

M 

X 

Sx 

Nx 

Tx 

0 

-0.24 

-0.66 

0 

0 

-0.25 

-0.66 

0 

0 

1 

0.23 

0.54 

44 

-12 

0.23 

0.54 

44 

-12 

2 

7.60 

18.00 

1720 

-307 

7.60 

18.00 

1700 

-350 

3 

5.60 

13.30 

1300 

-97 

5.50 

13.00 

1253 

-126 

4 

-4.50 

-10.00 

-1047 

0 

-5.00 

-12.00 

-1107 

0 

)tal 

8.60 

21.70 

1917 

-416 

7.98 

19.00 

1890 

-494 

Difference  expressed 

in  percentage  -7 

-12 

-1 

+18 

Note: 

Values 

of  Mx>Sx 

and  N 

A 

at  <J>  =  0 

and  N 

X(j) 

at  4= 

45° 
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TABLE  5.2 

VALUES  OF  STRESS  RESULTANTS 

AT  THE  SUPPORT 

M  at  <j>=0  N  at  <j>=0  N  ,  at  4>=1 5°  N,  at  4>=0° 

X  X  X(p  (p 

kips. ft/ft  kips  /ft  kips/ft  kips/ft 

Method  2 

-3.00  -118.2  -5.6  -29.95 

Method  3 

0  -118.2  -5.6  -29.95 

I  Q0H12M  X/V 
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FIGURE  5.7  CRITERION  FOR  MET  NOD  I  FIGURE  5.  8  CRITERION  FOR  METHOD 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

6.1  SUMMARY 

The  objectives  of  this  study  are: 

1.  to  obtain  usable  closed  form  expressions  for  stiffness  in¬ 
fluence  coefficients  for  closed  circular  cylinders  for  all 
harmonics  of  edge  loading  greater  than  zero  and, 

2.  to  establish  the  ranges  of  validity  of  these  expressions. 

Explicit  expressions  for  influence  coefficients  based 
on  the  exact  theory  are  very  unwieldy  and  their  evaluation  is 
tedious.  Therefore  approximate  solutions  discussed  in  lit¬ 
erature  were  considered.  Influence  coefficients  were  developed 
based  on  the  following  three  approximate  methods. 

Method  1  is  based  on  the  split  equations  of  Flugge 
for  low  values  of  m  and  the  solutions  obtained  are  termed  as 
the  'edge  effects'  and  the'basic  state.'  The  latter  is  re¬ 
placed  by  the  membrane  solution  within  certain  limitations.  The 
'edge  effects'  solution  considers  M  and  S  and  the  boundary 
conditions  correspond!" ng  to  w'  and  w  whereas  the  membrane  solution 
considers  N  and  T  and  the  boundary  conditions  correspondi ng 

X  X 

to  u  and  v.  The  stiffness  matrices  obtained  are  combined  to  get 
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the  matrix  [SM].  A  method  of  successive  correction  has  been 
proposed  to  consider  the  interaction  of  the  above  two  stress 
systems  and  obtain  the  final  stiffness  matrix  [SM].  Only 
one  cycle  of  correction  is  required  to  obtain  [SM]  from  [SM]. 

Method  2  is  based  on  Donnell's  theory  neglecting  the 
effects  of  the  far  edge  for  high  values  of  m. 

Method  3  is  based  on  the  split  type  stress  systems, 
the  generalised  plane  stress  and  the  plate  bending  for  very 
high  values  of  m. 

In  these  three  methods  explicit  expressions  were 
developed  for  all  elements  in  the  stiffness  matrices.  Criteria 
have  been  established  for  the  ranges  of  validity  for  these 
methods.  These  are  based  on  establishing  a  limiting  value  of 
m  corresponding  to  a  maximum  specified  tolerance  of  the  most 
critical  term  of  the  most  sensitive  stress  resultant  for  that 
method.  Examples  illustrating  the  application  of  these  stiff¬ 
ness  matrices  to  obtain  final  stress  resultants  are  given. 

6.2  CONCLUSIONS 

1.  The  solution  of  closed  cylindrical  shells  subjected 
to  non-axisymmetric  loading  may  be  conveniently  analysed  by  the 
method  of  stiffness  influence  coefficients.  Within  the  limit¬ 
ations  presented  it  is  possible  to  develop  simplified  expressions 
for  the  stiffness  influence  coefficients  in  closed  form  suitable 


for  practical  use. 
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2.  Expressions  for  the  roots  of  the  split  equations 
proposed  by  Flligge  for  the  characteristic  equation  for 
cylindrical  shells  can  be  simplified  by  equating  the  real  and 
imaginary  parts  (i.e.  x-|=r-|>  x2=y2 ^ roots  so  obtained 
are  in  close  agreement  with  the  exact  values  for  lew  values 
of  m. 

3.  Values  of  the  roots  expressed  in  terms  of  x-j  >  r-j  evaluated 
by  Donnell's  theory  agree  well  with  the  exact  ones  for  low  and  high 
of  m.  This  means  that  Donnell's  theory  also  can  be  used  to 
obtain  the  edge  effects  solution  in  Method  1. 

4.  For  low  values  of  m  the  roots  based  on  the  parameters 
X£>  are  sufficiently  small  to  be  considered  as  zero.  This 
results  in  permitting  the  replacement  of  the  basic  solution  by 
the  membrane  solution  in  Method  1.  The  maximum  value  of  m  is 
dependent  on  the  shell  geometry  and  increases  as  the  shell  becomes 
thinner  and  shorter. 

5.  For  high  values  of  m,the  stress  resultants  can  be 
evaluated  using  Donnell's  theory  and  neglecting  the  effect  of 
the  far  edge.  (Method  2).  The  minimum  value  of  m  also  depends 
on  the  shell  geometry  and  decreases  as  the  shell  becomes  thicker 
and  longer. 

6.  Based  on  a  maximum  specified  error  in  the  final  stress 
resultant  criteria  can  be  developed  to  establish  values  of  m 


for  any  given  shell  geometry  to  define  the  range  of  validity  of 
Methods  1  and  2. 
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7.  For  most  shells  of  dimensions  encountered  in  practice  within 
reasonable  limits  of  the  maximum  specified  error  the  criteria 

will  uniquely  specify  the  method  for  evaluating  the  stiffness 
influence  coefficients.  However  if  extremely  small  values  are 
specified  for  -the  maximum  error  it  is  possible  to  obtain  a 
few  values  of  m  which  are  between  the  limits  of  validity  established 
for  Methods  1  and  2.  For  these  values  of  m  Method  1  should  be 
used,  since  this  gives  conservative  values  of  the  stress  resultants 
for  these  harmonics. 

8.  In  the  analysis  of  cylindrical  shells  on  isolated 
supports  involving  very  high  values  of  m,the  membrane  stress 
resultants  converge  very  slowly  and  the  bending  stress  resultants 
converge  very  quickly.  The  solution  for  the  membrane  stress 
resultants  which  involves  a  large  number  of  terms  of  the  series 
can  be  simplified  by  using  the  split  solutions  of  Method  3. 
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APPENDIX  A 


A1 


A2 


A .  1  HQLAND1  S  CHARACTERISTIC  EQUATION 

For  a  complete  presentation  of  the  theory  refer 
refer  to  Hoi  and  ( 13) . 

The  eighth  order  characteristic  equation  according  to 
Hoi  and  is 


A^  -  4m J m2 -1  +  [1_  +  6m2(m2-l)]  A4 

-4m2(m2-l)]  m J m2-l  A2  +  m4  (m2-!)2  =  o 


A.  1 


where  k'  =  h‘ 


1 2(1 “U2)r2 


The  roots  for  this  equation  are  expressed  as 

4 


A 


1,2, 5, 6 


=  1  f2  A!'  [1  +  /  /l  +  e2  +  e  ±  1  (1  +Jj  l+£2  -e  )] 


>3  4  7  8  =±]n  4p  [i  -  In  +  z2  +  e  ±  i  (1  -//l+(?  -e  )] 


where ,  e  =  4  k 1  m  /  m  -  1 . 


A. 2 


A. 2  DONNELL'S  CHARACTERISTIC  EQUATION 


Donnell's  equation  is  given  in  the  form, 

4 

,  ?  2\  1  A* 

(A  -  m  )  +  p  =o 


A. 3 


This  equation  can  be  obtained  from  Hoi  and1 s  equation  by  making 
2  2 

m  -  1  =  m  .  The  roots  are  expressed  as, 


i  Ji  i 


X1 ,2,5,6  =  ±  '/2  I  Ik'  [U//]  +  +  eo  1  1(1  t/J1  +  eo  -  £o)]  ’ 


A3 


''“3,4,7 ,8 


where  e  =  4 
o 


m 


(A. 4) 


A .  3  SPLIT  EQUATIONS  OF  FLUGGE 

The  eighth  order  characteristic  equation  derived  by  Flugge 
has  been  split  into  two  fourth  order  character!- sti c  equations  in 
the  following  form. 


X4  -  2(2m2  -  u)x2  +  [1  -  u2  +  6m2 (rn2  -  1)]  =  o  • 

k 

[1  -  u2  +  6m2 (m2  -  l)]x4  -2m2  [2m4  -  (4  -  u)m2  +  (2  -  u)]x2  =  o*(A-5a,b) 

k 


The  large  roots  X-j  ^  5  5  are  given  by  the  first  equation  and  the 
small  roots  x0  „  0  by  the  second  equation. 

0  5  Tj/  >0 

EQUATION  (A-5a)  is  of  the  form, 

4  2 

X  +  bx  +  c  =  0  ,  where  > 

b=  -2(2m2  -  o)  and  c  =  1  -  u2  +  6m2(m2  -  1)  * 

k 

A2  =  [-|±L1/2'\/b2  -  4c]  • 


2  2 

As  4c  >  b  ,  (b  -  4c)  is  negative  > 
9 

X  =  [q -j  ±  i q 2 1  >  where  , 


A4 


^1>2  =  ±  (x-|  -  iun )  • 

x5,6  =  1  (x,  +  0,)  ’ 


Where,  x,  =  -  [  jq?  +  q/  +  q,]  /2  , 

_  1  r  m  2"  J/2  . 

pl  -  -  +  ^2  -  ^ 

c  >>  ^-  as  k  is  small  and  m  is  low.  Therefore  ignoring  q 
I  9  9* 

pared  to  +  q2  ? 

Xi  =rV21[1  -  u2  +  6m2(m2  -  1)]1/4  * 

1  k 

vi,  =  \}  /  2]  [1  -  u2  +  6m2  (m2  -  1)]  ^  • 

1  k 

It  is  seen  that, 


In  a  similar  manner  we  get  > 

x2  =[V2] 

y2=[V2] 


f  4,  2  ,>2 

m  ( m  -  1 ) 

(1  -  u2)  +  6m2 (m2  -  1) 
k 

[4/2  , n2 

m  (m  -  1 ) _ 

(1  -  u2)  +  6m2 (m2  -  1 ) 


1  com- 


(A-6a,b) 


(A- 7 ) 


(A-8a,b) 


A5 


It  is  seen  that, 

*2  =  y2  ' 

A. 4  EVALUATION  OF  a.  AND  3.  USING  THE  EXACT  THEORY 

- j- . -J - 

The  derivation  of  analytical  expressions  for  a.  and 

3.  using  EQUATIONS( 25-a ,b) ,  p.  227.  ref. 14  is  very  tedious. 

J 

For  particular  cases  they  can  be  determined  numerically.  The 
procedure  is  outlined  below. 

The  above  equations  can  be  reduced  to  the  form, 

(x2  +  s -j ) A  +  s2XB  +  (s3X  -  kX2)C  =  0  * 

s2AA  +  (s4X2  +  m2)B  +  (m  +  s5X2)  C  =  0  » 

where  the  following  substitution  have  been  made. 

sn  =  -  (1  -  u)m2  (1  +  k) 

1  2 

s9  =  (1  +  u)m 

L  2 

s9  =  u  -  k(l  -  u)m2 

j  2 

s-  =  -  (1  -  u)  (1  +  3k) 

4  2 

sc  =  -  (3  -  u)  km 

5  -y- 


>.-9a  ,b) 


A6 


Taking  any  root  A.  in  the  form  a  -  ib,  and  substituting  in  these 

J 

equati ons ,  we  get  > 

(x-j  +  ix2)A  +  (>'3  +  ix4)B  +  (xg  +  ix6)C  =  0  ? 

,(Xg  +  ix^)A  +  [x-j  +  iXg)B  +  (xg  +  ix-^)  C  =  0  (A-10a,b) 


where , 

xi 

a2  -  b2  ♦  s, 

11 

CNJ 

X 

2albl 

x3 

S3al 

X 

-p=> 

11 

S2bl 

x5 

s 3a-j  -  k[a-j(a-|2  -  b-,2)  -  2a-, b-j 2] 

x6 

s 3b -  k [2a2b ^  +  (a-,2  -  b-,2)b,] 

X 

II 

s4(a|2  "  +  ^ 

II 

CO 

X 

4a 1 b 1 

X 

CO 

II 

m  +  s5(a,2  -  b-, 2) 

II 

0 

X 

2s5albl 

EQUATIONS  (A-10a,b)  can  further  be  reduced  to  the  form  , 
A  +  (z-j  +  i  z  ^ )  B  +  (z3  +  iz4)C  =  0  » 


A  +  (z5  +  i z6 ) B  +  (z7  +  iz8)C  =  0 


9 


(A-lla,b) 


where, 


(x-,2  +  x/) 


Z1 

(x-|X3  -hx2x4) 

(x-, 2  +  x22) 

z5 

(X3X7  "  X4X8) 

!  2  ,  2x 

(x3  +  x4  ) 

z2 

= 

- ( Xp x3  -  x] x4 ) 

Z6 

= 

" (X4X7  ”  X3X8^ 

(X,  +  x2  ) 

(x3  +  x4  ) 

Z3 

= 

(xlx5  +X2X6> 

Z7 

= 

(x3xg  +  x4x10) 

(x]2  +  x22) 

(x3  +  x4  ) 

Z4 

= 

- (X2X5  "  X1X6^ 

Z8 

- 

"(X4X9  “  X3X 1 0 

f  2  ,  2x 

(x3  +  x4  ) 


as , 


Solving  EQUATIONS  (A-11a,b)  we  can  express  A  and  B  in  terms 


A  =  a  .C  and 
J 

B  =  3-C 

J 


where  a.  and  3.  are  as  follows, 
J  J 


a. 


( z  -|  +  i  z2 )  ( c-j  +ic2)  +  (z3  +  iz4)(c3  +  ic4) 

(c3  +  ic4) 


=  (Cj  +  ic2) 

(c3  +  ic4) 


The  following  substitution  have  been  made  in  the  above. 


A8 


C2 
C3 
C4 
c5 
C6 

a.  and  3.  are  complex  numbers.  The  real  and  imaginary  parts  of 

J  J 

these  can  be  written  as  follows, 


RefotjJ  (c^c^  +  CgC^) 


(c32  +  c42) 

Im[a  •  ] 

vj 

“(C4C5  -  CgCg 

,  2  x  2s 

(c3  +  c4  ) 

RqCBj] 

(C1 c3  +  C2C4) 

,  2  ,  2s 

(c3  +  c4  ) 

Im[ 3 j ] 

-(c1c4  -  C2C3 

,  2  x  2^ 
(c3  +  c4  ) 


From  these  quantities  values  of  a-j ,  a^,  a^,  a4  and  3^,  3^,  3^, 

3/.  in  terms  of  which  a.  and  3..  are  expressed  for  all  the  eight 
^  J  j 

roots  can  be  found  (see  p.  228,  ref.  14). 


z8  '  z4 


Z1  “  z5 


z2  -  z6 


■Clzl  +  C2Z2  "  C3Z3  +  C4Z4 


C1 z2  '  C2Z1  "  C3Z4  "  C4Z3 
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A . 5  EVALUATION  OF  a.,  3  .  FOR  THE  ROOTS  EXPRESSED  IN  TERMS  OF 

3  J 

X] jlXL]  for  the  split  equations  of  flugge 

The  following  assumptions  have  been  made. 


1.  X]  =  p-|  (see  SECTION  A. 3) 

2.  k  <<  1  and 

3.  1  «4x-|2  (see  SECTION  4.4.1  ) 


Analytical  expressions  for  the  root  x-j  -  iy-p  derived  using  the 
above  assumptions  and  the  procedure  outlined  in  SECTION  A. 4  have 
been  given  below. 


a.,-. 


c 


3 


1 


3 


2 


[-  u  +  (m2  +  4kx-r4)  ] 

0  2  2 ' 

2X-|  m 


2xl 


-[u  +  (m2  +  4k*-,4)  ] 


2  2 

2  x-j  m 


2xl 


[(2  +  u)  +  4kx-|4  ] 


m 


3 


-  [  (2  +  u)  x-| 


m 


]  _m 
2x 


A . 6  EVALUATION  OF  a  . .  USING  DONNELL 1  s  THEORY 

J  J 

For  a  complete  description  of  the  method  of  derivation  the 


A10 


reader  is  referred  to  Hoff  (12).  Only  the  final  expressions 
are  given  here. 


1 

1  2k 


an  = 


,2 


uy i  -  m  y-j 


xi  +  o 


2 


a0  -  -tti  1 2 


2  -  2k '  UX1  +  m  X1 


X1  +  0 


$1  -  -x-j  y  i 


wz(  2  +  2x2 

k  ( x i  +  y-j  ) 


32 


m  ^  [2  +  u  -  m2(Xl2  -  y-, 2)  ] 

2k ' 2  2  ,  2.2 

(Xi  +  y-j  ) 


a3 


1  2 

2k"1 2  t  uy2  "  m 
2  .  2 


X1 


+  y* 


a4  = 


2k 


,2 


ux2  +  m  x2 


2  +  2 
x9  +  y9 


e3  = 


-  x2y2  m 

.  ,2/  2  x  2) 

k  (xo  +  y9 


34  = 


m 

2k 


2,  2 


~2  [  2  +  o  -  m  (x2  -  y2") j 


(X22  +  P2Z) 


All 


A . 7  MEMBRANE  SOLUTION  FOR  SURFACE  LOADING 


For  a  complete  presentation  of  the  theory  proposed  by 
Flligge  refer  to  pp  107  to  132,  chapter  3,  ref  14.  In  this  ref¬ 
erence  Fliigge  has  given  EQUATIONS  (4)  for  the  stress  resultants 
and  EQUATIONS  (22)  for  the  displacements  in  the  dimensional 
form  using  the  sign  convention  shown  in  FIGURE  3.1.  The  co¬ 
ordinates  are  measured  from  the  centre  of  the  shell  as  shown 
in  FIGURE  A.l.  Taking  the  surface  load  in  the  form  pm  =  qm  cosm<j>, 

p  =p  =0  and  making  the  substitution  y  =  x_,  p=h,  and  Q=£_  the 
4  x  r  r  r 

above  equations  have  been  reduced  to  the  following  form  in  terms 
of  the  non  dimensional  parameters  y  and 


N 


EP  u 

EP  v 
EP  w 

EP  W 


Til 


m  q  y 
m 

TL(Q2-4y2)q 
8 


m 


m2(3Q2-4y2)  -  uy 
24 


m 


-  ( Q  2  -  4y  2 ) 
(Q2-4y2) 


(5Q2  :4y2)m2  +  2  +  u 
.  384  8 


m 


(5Q2-4y2)m2  +  ]/4  m2  +  1 
l  384 


in 


-8y 


r,cr2  n  2s  2,1,.  ,  fr2  .  2s  2 

(5Q  -4y  )  m  +  /4  +  (Q  -4y  )  m 

384  384 


2 

m  q 


m 


To  get  the  absolute  values  of  the  stress  resultants  N  ,  N^  and  Nx 
and  the  displacements  u" ,  V  and  w,  these  expressions  are  to  be 
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multiplied  by  the  radius  r. 

A . 8  FOURIER  SOLUTION  FOR  EDGE  LOADING  BY  MEMBRANE  THEORY 

For  a  detailed  description  of  the  method  refer  to 
section  3. 2. 2. 3,  chapter  3,  ref  14.  In  this  reference  Flugge 
has  given  the  homogeneous  solutions  for  the  displacements  given 
by  EQUATIONS  (25)  and  (26)  in  the  dimensional  form.  Making  the 
same  substitutions  as  explained  in  SECTION  A. 7  these  equations 
can  be  reduced  to  the  following  non  dimensional  form, 


EP  u 

= 

A  y  + 
2 

By  +  c 

EP  v 

A[m* 2y3 

-  2(l+u)my]  +  Bmy2  +  Cmy  +  D 

6 

2 

EP  w 

- 

-A[m2y3 

-  (2+u)my]  -  B[m2y2  +  u]  -Cm2y 

6 

2 

EP  w' 

- 

nr  2  2 

-A[m  y 

-  (2+o)m]  -  Bm2y  -  Cm2 

2 

To  get  the  absolute  values  of  the  displacements  u,  v  and  w  these 
expressions  are  to  be  multiplied  by  r. 
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£*>  CONST 


A/x  +  dx 


2x 


d0 


N$  *  dilt  cJd 

d& 
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